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ABSTRACT 

We identify a "supernova rate problem": the measured cosmic core-collapse supernova rate is a factor of 
^ 2 smaller (with significance ^ 2a) than that predicted from the measured cosmic massive-star formation 
rate. The comparison is critical for topics from galaxy evolution and enrichment to the abundance of neutron 
stars and black holes. We systematically explore possible resolutions. The accuracy and precision of the star 
formation rate data and conversion to the supernova rate are well supported, and proposed changes would 
have far-reaching consequences. The dominant effect is likely that many supernovae are missed because they 
are either optically dim (low-luminosity) or dark, whether intrinsically or due to obscuration. We investigate 
supernovae too dim to have been discovered in cosmic surveys by a detailed study of all supernova discoveries 
in the local volume. If possible supernova impostors are included, then dim supernovae are common enough 
by fraction to solve the supernova rate problem. If they are not included, then the rate of dark core collapses is 
likely substantial. Other alternatives are that there are surprising changes in our understanding of star formation 
or supernova rates, including that supernovae form differently in small galaxies than in normal galaxies. These 
possibilities can be distinguished by upcoming supernova surveys, star formation measurements, searches for 
disappearing massive stars, and measurements of supernova neutrinos. 

Subject headings: galaxies: evolution - galaxies: starburst - stars: formation - supernovae: general 



1. INTRODUCTION 

Core-collapse supernovae (CC SNe) are extremely im 
portant to many areas of astrophysics. They are re 
sponsible for the major ity of heavy elements (e.g 



iMatteucci & Greggiolll986 ), are a ssociated with dust produc- 
tion (e.g., i Todini & Ferraral|200 l'). and could dominate winds 
and feedback in galaxy formation (e.g., iFerrar a & Ricott^ 
l2006h . Observations of CC SNe and their progenitors test our 
understanding of stellar evolutio n (see, e.g., the progenitor- 
SN map of Gal - Yam et al. 1 I2007I) . and the extreme densities 
and temperatures reached in CC SNe offer the opportunity 
to study the physics of weakly interacting particles, forma- 
tion of compact objects, and related nuclear physics (see, 
e.g., Raffelt 1990, 2000: Kotake et alJ200d : ljanka et alJ2007l: 
lEattimer & Prakash.2007i, for recent reviews). However, there 
is much that is not yet understood. 

One of the outstanding questions is the initial conditions 
correspo nding to optically luminous, dim, and dark CC SNe 
(see, e.g.. |Janka et alj|2"007l) . It is expected that some CC SNe 
will be intrinsically dim or even dark in the optical bands. 
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One possibility is the collapse of stars near the SMf;^ t hreshold 
via an electron capture trigg er (iMivaii et al.lll980t iNomotol 
Il984l: IPoelarends et al.ll2008l) . Another possibility is the col- 
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Fig. 1 . — Comoving SNR (all types of luminous core collapses including 
Type II and Type Ibc) as a function of reds hift. The SNR predicted f rom 
the cosmic SFR fit and its supporting data iHopkins & Beaconill2006l) . as 
well as that predicted from the mean of the local SFR measurements, are 
plotted and labeled. The fit to the measured cosmic SNR, with a fixed slope 
of (1 + z)"^'^ taken from the cosmic SFR, is shown with the uncertainty 
band from the LOSS measurement. The predicted and measured cosmic SNR 
are consistently discrepant by a factor ~ 2: the supernova rate problem. 
However, rates from SN catalogs in the very local volume do not show such 
a large discrepancy (see Figurejs]. 
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lapse proceeding directly to a black hole. Theoretical stud- 
ies indicate that p rogenitors w ith M > 40 Mq may promptly 
form black holes ("Frved llQQQl) . with additional dependencies 
on the stellar meta llicity and rotation of the progenitor core 
(iHegeret al.ll2003h . CC SNe that are intrinsically luminous 
can also appear as optically dim CC SNe if they are heavily 
enshrouded by dust. Due to their darkness, these collapses are 
difficult to study directly. 

Precise measurements of the cosmic star formation rate 
(SFR) and the cosmic CC SN rate (SNR; we include all Type 
II and Type Ibc SNe) provide new information about CC SNe. 
The SFR and SNR encode the birth rate of stars and the death 
rate of massive stars, respectively. Since the massive stars 
that give rise to CC SNe have cosmologically short lifetimes 
^ 30 (A//8A/0)"2-^ Myr, the cosmic SNR is expected to 
follow the same evolutionary trend in redshift as the cosmic 
SFR. Second, the SFR data show that massive stars are clearly 
present, so they must be dying even if not all of these deaths 
are optically luminous. Because the measured SNR are sensi- 
tive to only optically luminous CC SNe, the normalization of 
the SNR with respect to the SFR provides information on the 
frequency of massive stars' fates as either optically luminous 
or dim (and dai'k) CC SN^l 

Early attempts to compare the evolution and normaliza- 
tion of the cosmic SFR and cosmic SNR were in c onclu- 
sive due to their large uncertainties. IPahlen et al ] (l2004 
found reasonable agreement between their measured SNR 
at z fti 0.3 and z sa 0.7 and a compilation of SFR data 
by iGiavalisco et al.l (l2004l) . Us ing an improved compila- 
tion of modern cosmic SFR data, Eo^ ins & BeacomI (120061) 
noted that the measured cosmic SNR (Cappe llaro et al.ll 19991 
iDahlen eTal. .2Q04; Cappellaro et al. 2005) were somewhat 
lower than those predicted from the SFR . Sim ilar conclusions 
were reached by iMannucci et al] (l2007h and iBotticella et al.l 
(12008.) . 

In recent years, measurements of the cosmic SFR and 
cosmic SNR have rapidly improved. The cosmic SFR has 
been measured using multiple indicators by many compet- 
ing groups. The accuracy and precision of the cosmic 



SFR has been documented (e.g., iHopk ins & BeacomI 



and is supported by recent data (e.g .. iPascale et al.l 



2006) 



2009 



iRuiopakarn et al. I I2OIOI; iLv et al.ll201 lllBothwell et aliboi lb 
The Lick Observatory Supernova Search (LOSS) has recently 
published the best measurement of the cosmic SNR at low 
redshifts, using CC SNe collected over many years of system- 
atical l y surveying gala xi es within ^ 200 M pc CLeaman et alj 
l201lt Li et al.l I20 11 alibi: iMaoz et all 1201 ll) . The Supernova 
Legacy Survey (SNLS) has published the most precise SNR 
measurement at higher redshifts, using a large sample of CC 
SNe c ollected in their e xtensive rolling search of four deep 
fields teazin et alj |2009'). 

Based on the latest data, it has become clear that the mea- 
sured cosmic SFR and the measured cosmic SNR both in- 
crease by approximately an order of magnitude between red- 



' Typically, dim (low-luminosity) CC SNe are loosely defined as CC SNe 
falling at the lowest end of the CC SNe luminosity function; here, we define 
them as dimmer than M 
results; see Section|3.6| 



-15 mag. The two definitions yield similar 



shift and 1, confirming our expectation that t he progeni- 
tors of CC SNe are sh ort-lived massive stars (e.g.. lBazin et alj 
l2009HLi etaDl2011bl) . On the other hand, the comparison of 
the normalizations of the latest SFR and SNR data has been 
left for future work. We perform this here for the first time. 
As illustrated in Figure [T] the SNR predicted from the cosmic 
SFR is a factor of ^ 2 larger than the cosmic SNR measured 
by SN surveys; we term this normalization discrepancy the 
"supernova rate problem." Both the predicted and measured 
SNR are of optically luminous CC SNe, so the two can be 
directly compared. The lines in Figure [Hare fits to the SFR 
and SNR data, respectiveljH The discrepancy persists over 
all redshifts where SNR measurements are available0. 

The nominal uncertainties on the fits (shaded bands) are 
smaller than the normalization discrepancy, and the signifi- 
cance of the discrepancy is at the ^ 2a level. At high red- 
shift, where the uncertainties of the SNR measurements are 
largest, the statistical significance is weaker However, it is re- 
markable how well the cosmic SNR measurements adhere to 
the expected cosmic trend — much better than their uncertain- 
ties wo uld suggest. Indeed, the measurements of Dahlen et al] 
(l2004 have been supporte d by recent unpublis hed results and 
with reduced uncertainties (iDahlen et al.ll20l6l) . We therefore 
consider the fits to be a good representation, i.e., the super- 
nova rate problem persists over a wide redshift range. We sys - 
tematically examine resolutions to the supernova rate prob- 
lem, exploring whether the cosmic SNR predicted from the 
cosmic SFR is too large, or whether the measurements under- 
estimate the true cosmic SNR, or a combination of both. 

In Section 2, we describe the predicted and measured cos- 
mic SNR in detail and substantiate the discrepancy. In Section 
3, we discuss possible causes. In Section 4, we discuss our re- 
sults and cautions. We summarize and discuss implications in 
Section 5. Throughout, we adopt the standard ACDM cos- 
mology with = 0.3, Q.A — 0.7, and Hq — 73 km s^^ 
Mpc"^ 

2. NORMALIZATION OF THE COSMIC SNR 

The cosmic SNR is calculated from the cosmic SFR using 
knowledge of the efficiency of forming CC SNe. The most 
recent SFR is traced by the most massive stars that have the 
shortest lifetimes. The primary indicators of massive stars — 
Ha, UV, FIR, and radio — are routinely used, with dust cor- 
rections where necessary, to study the populations of massive 
stars. However, since the total SFR is dominated by stars with 
smaller masses, the SFR derived from massive stars must be 
scaled upward according to the initial mass function (IMF); 
for example, for a given massive stellar population, an IMF 
that is more steeply falling with mass will yield a larger to- 
tal SFR compared to a shallower IMF. The scaling is done 
with the use of calibration factors derived from stellar popu- 
lation synthesis codes that calculate the radiative output from 

Technically, the SNR line shown is not a fit, but is a conservative estimate 
based on the SNR measurement of LOSS; see Section|2] 

' However, in the local < 25 Mpc volume, the SNR derived from SN 
catalogs does not show such a large discrepancy, sup porting earlier claims 
that the true c osmic SNR is as large as predicted fe.g., lHoriuchi et ai]l2009l ; 
lBeaconj20"Toh . 
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TABLE 1 

Summary of Cosmic SNR Measurements 



Redshift 


Rate (10"'' yr-i Mpc'^) 




Cadence (days) 


JVcCSN 


Host Extinction 




(mag) 




(mag) 


Reference 
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0.21 


-, 0-1+0.49 1^+0.48 N 
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~ 23 
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-16.7 


3 


0.26 


-I O-I+0.91 
-^■°-'--0.79 


11300 




19.5 


N 


mv 


~ 23 


Mv ~ 


-17.1 


4" 


0.3 


1 nt:+0.34/+0.34-, 
J--OO-0. 341-0. 34^ 


Volumetric 


~ 7 


117 
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rrii 


~ 24 
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-16.4 


5 


0.1-0.5 


9 oc + 1.00|'+0.85n 
^■f'a-0. 85^.-2. 11^ 


Volumetric 


~ 45 


6 


Y 


rriz 


~ 26 


Mr ~ 


-14.8 


6 


0.5-0.9 


^■^^-1.20 '--2.95-' 


Volumetric 


~ 45 


10 


Y 


mz 


~ 26 


Mv ~ 


-16.0 
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References. ^ (l) ICappellaro et ai] Jl999ft : (2) lLi et al.l f201 Ibl) ; (3) lBotticella et al.N2008l) : (4) ICappellaro et al.l | |2005| ): (5) lBazin et all <2009l) : (6) 
lDahlen"et '^ f2004l) . 

Note. — Systematic error estimates are given in parentheses. The limiting absolute magnitudes are calculated from the limiting apparent magnitude 
at the distance of the data point. The scope of the rate measurements depend on a combination of limiting magnitude, SN luminosity function, and dust 
connection (see the text). 

" We adopt the rate per B-band luminosity to volumetric rate conversion used by Bazin et al. (2009). 
'' Includes the CC SN sample of Cappellaro et al. (2005). 



a population of stars following an assumed IMF (see, e.g., 
iKennicuti 19981) . 

We adopt the dust-corrected SFR compilation of 
iHopkins & BeacomI (,2006). Their dat a are well fit by a 
smoothed broken power law of the form (lYiiksel et al.ll2008h 



to Mmax, resulting in a SNR of 



RsNiz) = f>*{z) 



(2) 



P*iz) = Po 



(1 + ^)"" + 



1 + z 
B 



b7] 



1 + z 
C 



iX) 



where B = (1 + (7 = (1 + + 

zif-^l". We adopt = 0.016 /173 Mq Mpc^^ yr-i for 
the cosmic SFR at z = 0, as well as the parametrization 
a = 3.4, h = -0.3, c = -3.5, zj = 1, Z2 = 4, and 
r\ = —10. These choices are applicable for the Salpeter A 
IMF, whic h is a modified Salpeter IM F with a turnover be- 
low IMq dBaldrv & Glazebrookll2003l) . The scaling from a 
Salpeter IMF is « 0.77. The Icr uncertain ty on po is approx- 
imately ±25% dHopkins & Beacomll2006l) . The fit is in good 
agreement with a range of recent SFR measurements over a 
range of redshifts, for example using Ha (Ly et al. 20 1 j), UV 
jSalim et al.ll2007b. IR (jPascal e et al. 20091: iRuiopakarn et al.1 
I2OIOI) . and X-ray (IWatson et al.i .2009). Most recently, a thor- 
ou gh study of the local S FR at z w 0.05 has been performed 
bv lBothwell et alj(l201lh using combined observations in UV 
and IR. Adjusting to our chosen cosmology and IMF, the de- 
rived SFR is 0.0193 ± 0.0012 /173 Mq Mpe'^ yr-\ in good 
agreement with our parametrization. At ±6%, the uncertainty 
is a great improvement over many previous measurements and 
our adopted uncertainty on po- 

The comoving volumetric SNR is determined by multiply- 
ing Equation ([T]l by the efficiency of forming CC SNe. This 
is the number of stars that eventually explode as CC SNe per 
unit stellar mass formed in a burst of star formation. It is 
largely governed by the mass range for CC SNe, from A/min 



where ijj{M) is the IMF, defined over the main-sequence mass 
range 0.1-100 Mq. The IMF is defined such that ^{M)dM 
gives the number of stars in the mass range A/ to M + dM. 
Due to the steeply falling nature of the IMF, the lower mass 
limit Mniin is the most important parameter Note that this 
selection of the relevant mass range is effectively the inverse 
process of the scaling from the massive-star SFR to the total 
SFR. In fact, the stellar mass range probed by the SFR indi- 
cators is comparable to the mass range giving rise to CC SNe. 
Thus, variations in the IMF should have only a small effect on 
the predicted SNR. 

The predicted cosmic SNR is shown in Figure [T] We as- 
sumed canonical parameters for optically luminous CC SNe, 
Af,„i„ = 8Mq and Af,„ax = 4OM0. The SFR to SNR conver- 
sion coefficient is then O.OO88/M0, yielding Rs^iz = 0) w 
1.4 X lO"** yr~^ Mpc^^. The uncertainty band shown is the 
Icr uncertainty in the SFR fit of iHopkins & BeacomI (|2006|) . 
The actual data of the compilation is also shown, similarly 
converted to a SNR: th e dust-corrected UV from Sloan Digi- 
tal Sky Survey (SDSS: Ba ldrv et al.ll2005b. Galaxy Evolution 
Explo rer (GALEX: Arnouts et al.ll2005[ ISchiminovich et al.l 
l2005l) . C 0MB017 iwolf et alJl2003l) . and Ha-derived mea- 
surement (' Hanish et a l. 20061 The mean of the local SFR 
measurements (Section |372] i is also converted and shown. 

There are two main approaches for collecting CC SNe and 
measuring the cosmic SNR. In the first, the same patch of 
sky is periodically observed, locating CC SNe within a vol- 
ume, limited only by flux. In the second, a pre-selected sam- 
ple of galaxies is periodically observed, and the rate per unit 
galaxy size (mass or light) is converted to a volumetric rate 
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using the galaxy mass or luminosity density. In the former, 
the completeness of the CC SN sample is readily definable 
and the volumetric rate is derived directly, but it requires a 
sufficiently wide or deep field to collect CC SN statistics. In 
the latter, there is a bias against CC SNe in small galaxies, 
but, for local surveys, it maximizes the CC SN discovery rate 
given the closer target volume. In Table [T] we summarize 
measurements of the cosmic SNR in the literature, showing 
the rate, the number of galaxies sampled (where appropriate), 
the cadence, the number of CC SNe used for the rate analysis, 
whether host galaxy extinction corrections are made, and the 
limiting magnitudes of the surveys; the survey characteristics 
vary considerably. 

The limiting magnitude of SN surveys is typically defined 
as th e magnitude at wh i ch the detection efficiency is 50% 



(e.g .. IDahlen et alJ 



For 'Botticella et al 



2004t iBazin et al ] 120091: iLi etatibOIlbh . 
(l2008l)T we quote the limiting magnitude 



for a seeing of 1 arcsec, their average. The limiting abso- 
lute magnitude is estimated from the limiting a pparent mag- 
nitude and the distance of the data point (for ^ Dahlen et al.l 
12004, we use the redshift that divides the distance bin in two 
equal volumes). This estimate is only indicative, because dim- 
mer CC SNe can be discovered at smaller distances. Fur- 
thermore, SNR measurements adopt a CC SNe luminosity 
function, usually but not always derived from the SN survey 
data itself, to correct for missing dimmer CC SNe. The SNR 
measurements should therefore be treated as measurements 
of CC SNe within the luminosity limits of the SN luminos- 
ity function adopted. We show apparent magnitudes in the 
observed-frame band and the absolute magnitudes in the su- 
pernova rest-frame band. For LOSS, we quote the limiting 
magnitude of their volume-limited subsample and not their 
full sample, since the former is used to correct the latter 

The most r e hable SNR measuremen t is that by LOSS 
jLeaman et all l201 ll : [Uet al. 201 la b; iMaoz et alj boi ih : 
they have excellent cadence and the largest published SN sam- 
ple. They use a large subsample of 101 volume-limited CC 
SNe to construct a modem SN luminosity function. It is com- 
plete to Mr ^ — 15 mag, although it is co rrected to go dim- 
mer (the LOSS-LF sample: lLi et aPbOl lal) . Their uncertainty 
of approximately ±30% is dominated by systematics, and 
originates mainly from limited knowledge of the z « galaxy 
luminosity density needed to deri ve the volumetric rates from 
their rate per galaxy luminosity dLi et al.ll201 Ibb . Not sur- 
prisingly, this is comparable to t he uncertainty o n the cos- 
mic SFR normalization at z w dHopkins & Beacom 200 6). 
At higher redshifts, the most precise measurement is that of 
SNLS, which is a volume monitoring search with excellent 
cadence and CC SNe statistics (Bazinetal. 2009). The lu- 
minosity function of SNLS is in good agreement with that 
of LOSS-LF, although LOSS goes approxima tely 1 . 5 ma g 
fainter on the supernova luminosity function (lRichll2010l) . 
Their result has been corrected for intrinsic extinction due to 
host inclination. 

Measurements of the cosmic SNR are shown in Figure [T] 
Those with host galaxy extinction corrections are indicated 



by filled symbols dCappellaro et al. 



iBotticellaetai] 120081; iBazin et al. 



1991 IDahlen etaDl2004t 



20091) . while those with- 



out are indicate d by empty symbols (ICappellaro et al.ll2005l 
iLi etalJl20TTbl) . The error bars combine the statistical and 
syste matic uncertai nties in quadrature. Unpublished results 
from D ahlen et al. OOlO) support their previous measure- 
ments (Da hlen et al.i r2004). The high redshift measurements 
are further supported by recent results from the Subaru Deep 
Field whi ch lie very close a lthough with somewhat larger er- 
ror bars (iOraur et al.ll20I ll) . Preliminary results from SDSS 
on the SNR at z < 0.2 also lie on the trend of the other rate 
measurements (iTavlor et al.ll2009b . 

A fit to the SNR measurem ents, excluding the measure- 
ment by lCappellaro et all (Il999 ) and assuming a slope identi- 
cal to the cosmic SFR, yields a best-fit normalization 0.70 x 
10^** yr~^ Mpc""^ with Xmin ~ 0.5 for 5 degrees of free- 
dom. The normalization is necessarily low due to the fact that 
some data have not been corrected for host-extinction. The 
is somewhat low, reflecting the nature of the large system- 
atic uncertainties. The \<j confidence region corresponding to 
< + Ay ^, where Ax^ = 1-00 for a 1 parameter fit 
(e.g..^ Avnilll976l) . is smaller than the systematic uncertainty 
of LOSS. To remain conservative, we show the combined sta- 
tistical and systematic uncertainty of LOSS and not of the fit. 
As a check, we also let the slope vary and find comparable 
results: a best-fit normalization of 0.72 x 10~^yr~^ Mpc""^ 
and a slope of 3.3, with Xmin ~ 0.5 for 4 degrees of freedom. 

The comparison of the predicted and measured cosmic 
SNR, and the relative sizes of the uncertainties, demonstrate 
two key points: they evolve similarly in redshift, and there is 
a systematic normalization mismatch. 

3. POSSIBLE EXPLANATIONS 

3.1. Is the Cosmic SFR too High? 

We address t he scatter in and check the n ormalization of the 
cosmic SFR of iHopkins & BeacomI (|2006|) . In particular, the 
scatter mainly originates from the different indicators used to 
derive the SFR, but also the intrinsic uncertainty in the cali- 
bration of each indicator The normalization is strongly influ- 
enced by the choice of IMF as well as dust corrections. Below, 
we discuss these in turn. 

The SFR compilation includes measurements from Ha, 
UV, FIR, and radio indicators. The scatter in the SFR fit 
is dominated by the scatter in the measurements from these 
different groups ; it is typi cally ±30% in the redshift range 
< z < 1 (Hop kins & Be acom 2006). At the lowest redshift, 
the scatter approaches the intrinsic uncertainties in the calibra- 
tions of each indicators. For the bulk of the galaxies the con- 
tinuous star formation approximation (i.e., star formation re- 
mains constant on time scales longer than the lifeti mes of the 
dominant UV-emitting massive stars, ^ 60 Myr; iKennicut^ 
Il998l) holds, so that the calibration uncertainties are ±10%- 
20%. If the assumption is incorrectly applied, as would be the 
case for young starburst galaxies, the SFR would be underes- 
timated by ^ 30%. Therefore, the scatter in the SFR derived 
from a single indicator and the scatter among different indica- 
tors are generally at the tens of percent level. 

The choice of IMF plays an important role in the normal- 
ization of the SFR, because the IMF is used to scale the ob- 
served massive-SFR to the total SFR. Adopting instead of the 
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Salpeter A IMF a classic Salpeter IMF or a flatter Baldry- 
Glazebrook IMF, the total SFR normalizati on varies at the 
±30% level (e.g.. 'Hopki ns & Beac om" 2006). However, it is 
important to stress that this does not affect the predicted SNR 
normalization anywhere near as strongly, because the massive 
stars that are used to make the SFR measurements are close 
in mass to the progenitors of CC SNe. In other words, an 
increase or decrease of the SFR due to steeper or shallower 
IMF is nearly completely canceled by a respective decrease 
or increase in the fraction of massive stars. The difference be- 
tween the predicted SNR using the Salpeter A and the B aldry- 
Glazebrook IMF is less than 10% (Horiuchi et al.ll2009l) . 

While extinction by dust raises the possibihty for the SFR 
to be over-corrected for dust, it is not realistically possible 
that the SFR in the most important redshift range (0 < z < 1) 
is over estimated by a factor of ^ 2. Since interstellar dust 
absorbs UV emission and re-emits in the FIR, the UV-derived 
SFR must be corrected. In our adopted SFR compilation, the 
UV-deri ved SFR are dus t -correc ted by adding the FIR-derived 
SFR of iLe Floc'h et alj (l2005h at their respective redshifts. 
The accuracy of this exti nction correction has been discussed 
using individual galaxies jlglesias-Paramo eT al. 2006) and b y 
comparisons with independ ent Ha-derived SFR (iBellll2003b : 
see also Hopkins ( 2004j) an d iHopk ms & BeacomI B006l) . The 
SFR corrected for dust in this way are consistent with recent 
SFR measurements in many wavebands and dust-correction 
methods (see Section |2|i. In the range z < 0.4, where the 
SNR normalization discrepancy is most significant, dust cor- 
rection typically increases the UV-derived SFR by a factor of 
^ 2. Thus, to explain the supernova rate problem by lowering 
the SFR, one would require almost no dust at these distances, 
which is inconsistent with our understanding of t he FIR uni- 
verse (iHauser & DweklbOOUlTagache et allboOSl) . 

At higher redshifts (0.5 < z < 1) the dust correction to 
the UV data become larger. By z « 1, more than 80% of 
the total SFR comes from the FIR. Eventually, the dust opac- 
ity in galaxies become high enough that the FIR luminosity 
itself measures the bolometric luminosity of the star forma- 
tion activity. To reduce the SFR by a factor of ~ 2 under 
these conditions requires one to either have severe contami- 
nation of the FIR emission from non-star-forming sources or 
that the FIR to SFR calibration is incorrect. While it is known 
that there is a "cirrus" component in the FIR that is associated 
with more extended dust heated by a population of older stars, 
their relative contribution varies substantially from galaxy to 
galaxy. For spirals it can be as high as 50%-70%, but for the 
star forming galaxies that d ominate at high redshifts it is ex- 
pected to be much smaller dLonsdale Persson & Heloull 19871 
iKewlev et all2002h . The calibration uncertainties are ±10%- 
20% as described above. 

Finally, the integrated SFR nor malization can be checked 
against observable quantities (Mad au et al.lll998l) . for exam- 
ple the stellar mass density and the extragalactic background 
light (EBL). The stellar mass density has been measured up to 
redshifts of a few and are typically a factor of 2 smaller than 
those pred icted from the SFR assumin g the Salpeter A IMF 
(see, e.g.. iHopkins & Beacomt .2006: Wilkins et alj [2008b). 
While a potential explanation of this discrepancy is to reduce 



the cosmic SFR normalization by a factor of ~ 2, i.e., at the 
same time also solving the supernova rate problem, this is un- 
likely to be the correct explanation. First, the stellar mass den- 
sity is dominated by low-mass stars rather than the high-mass 
stars constrained by the SFR, leaving the comparison to re- 
veal more about the IMF slope rather than the SFR normaliza- 
tion. For example, an IMF with a high-mass slope shallower 
(—2.15; equivalent to the Baldry-Glazebrook IMF) than the 
Salpeter slope (—2.35) reconciles the discrepancy at z < 0.4 
(Wilki ns et al. 2008a). Second, it has been argued that ob- 
servat ions could be underestimating the stellar ma ss density 
(e.g.. iNagamine et al.ll2004l:ISomerville et al.l200lb . Third, as 
we discuss below, such an explanation creates a tension in ex- 
plaining the EBL. 

The EBL is powered by moderately high mass stars (half of 
the stellar EBL is powered by stars with masses M > 3Mq), 
so that the IMF dependency is modest. Approximately 70% 
of the EBL arises from the z < 1 range we are most in- 
terested in, with only a few percent contribution from non- 
nucleosynthesis energy sources such as active galactic nu- 
clei (AGNs; Hopkins et al. 2006). This makes the EBL a 
useful probe of the integrated (0 < z < 1) SFR normal- 
ization. The minimum EBL has been derived indirectly by 
counting visible galaxies, whil e the maximum EBL ha s been 
directly measured (see, e.g^ jHauser & Dwek 2001]). We 
summarize the findings of IHoriuchi et al.. (2009i) here. In- 
tegrating from the FIR to UV bands yields a total EBL in 
the range 50-100 nWm~^ sr~^, comparable w ith oth er pre- 
viously evaluated estimates (e.g., Fardal et al] l2007h . The 
EBL has recently been constrained by the observations of 
gamma rays from di stant bla zars to be between the se ex- 



tremes jAharonian et al. 20 06: i Albert et al. 2008: Abdo et all 
l20IOt IOrretal.ll20Ill) . We adopt a nominal total EBL of 
73 nWm~^ sr~^ that respects the gamma-ray constraints. If 
secondary gamma rays produced en route by beamed cos- 
mic r ays are important (lEssev & Kusenkoll2010t lEssev et al] 
I20T1I) . the EBL could be even highe r. In comparison, 
the result calculated from the SFR of iHopk ins & BeacomI 
('2006"), using the PEGASE.2 p opulation synthesis code 



( 2UUb) , usmg tne j-i^UASfc.z p o| 
tFioc & Rocca-Volmerangd Il997l) . 



for flie Salpeter A IMF and 78^24 nWm'^ sr^i for the 
Baldry-Glazebrook IMF. Therefore, the predicted and mea- 
sured total EBLs do not allow a factor of ^ 2 reduction of the 
SFR. 

In summary, the uncertainties associated with SFR mea- 
surements are generally at the tens of percent level. The in- 
tegrated cosmic SFR normalization has been cross-checked 
with the EBL at a similar precision, and does not allow the 
cosmic SFR to be decreased enough to explain the supernova 
rate problem. 

3.2. Does the Local SFR Differ from the Cosmic SFR? 

The local (< 100 Mpc) SFR sets the birth rate of stars 
nearby and has important implications for our study of the 
local SNR. Here, we discuss whether the local SFR measure- 
ments are different from the cosmic SFR. Since the SFR are 
densities, they should not depend on distance, except for evo- 
lution with redshift, which is modest within 100 Mpc. We 
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address the validity of the comparison in Ught of expected 
cosmic variance and bias due to missing massive galaxies in 
small volumes, before discussing the results. 

Cosmic variance falls approxim ately as distance to the 
power — 2 (|Peacock&Doddslll994 and is minimal by sev- 
eral tens of Mpc. Normalizing the density powe r spectrum to 
the se ven-year WMAP result of as = 0.809 jjarosik et al.l 
I2OI lb . the density fluctuations within 30 Mpc should only 
be 0.1 of the mean density. Variations in smaller vol- 
umes are more likely. For example, it is not surprising that 
the galaxy S-band luminosity density in the 8 Mpc volume 
is a factor of 1.7-2.0 times the global luminosity density (as 
derived from the SPSS a nd the Millennium Galaxy Catalog; 
iKarachentsev et alj|2004 . 

Given the smaller volume, the rarest and most massive 
galaxies would be missing. However, the effect of missing 
their star formation is minimal. While massive galaxies have 
higher SFR, massive galaxies are more rare and have lower 
star formation per mass, so that the contribution to the total 
SFR peaks at galaxy mass ^ 6 x 10^" Mq. Consider the 
volume within 40 Mpc. Integrating the SDSS galaxy stel- 
lar mass function, the number of galaxies more massive than 
4 X 10" M0 is expected to be < 1. The fraction of the SFR 
misse d by excluding galaxies above 4 x 10^^ Mq is less than 
10% ( Brinchmann et al ] l2004l) . Therefore, the volume within 
several tens of Mpc is a fair representation of the global SFR. 

In Figure |2] we show SFR m easurements i n the 
local volume, derived f r om Ha ( Gallego et al. 1995 
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2008h. FIR dYun et all l200lh . radio dSerieant et al.l 12002 
Condon et al. '2002"), and combination of UV and IR 



jMartin et al. 2005; Bothwell et alj|201li) . The measurements 
are corrected for dust and have been scaled to our chosen cos- 
mology and IMF. At distances just beyond those shown, a 
recent thorough study using a combined UV and IR dataset 
find 0.0193 ± 0.0012 /173 Mq Mpc^^ yr^^ (at z « 0.05), in 
excellent agreement with the cosmic SFR and with small un- 
certainty (iBofliweUet al.ir2 011). 

The SFR within 11 Mpc has been measured by the 
11 Mpc Ha a nd U l traviolet Galaxy Sur vey (IIHUGS; 
iKennicutt etal] 120081; iBofliweU et al] 1201 ih . The value. 



0.0262 



+0.0044 
-0.0035 



/i73 Mq Mpc~'^ yr~ , shows an overdensity 
with respect to the cosmic SFR. The degree of overdensity 
is approximately a factor of 1.7 and is within expectations 
of cosmic variance. It is smaller than but similar to the en- 
hancement seen in the galaxy B-band luminosity density. The 
remaining measurements at larger distances are compatible 
within errors with the cosmic SFR. The mean of the mea- 
surements, omitting the IIHUGS measurement because of 
the known overdensity, is 0.017to oo5 Mq Mpc"^ yr"^ This 
is in good agreement with the cosmic SFR. The uncertainty 
range has been conservatively taken to be the range of the 
SFR measurements, because some of the measurements are 
in disagreement with each other, suggesting the true uncer- 
tainties are larger than reported for some individual measure- 
ments. We conclude that the local SFR data do not support 
large deviations from the cosmic SFR. 
In summary, SFR measurements in the local volume are in 
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Fig. 2. — Local SFR density as a function of distance. Measurements are 
shown at t he mean distances of their respecti ve galaxy samples , except for 
IIHUGS tBothwelletai]|20IIh and HaGS jj ames et aT]|2008l) which are 
measurements within fixed distances and are shown at their distance limits 
(horizontal arrows). The SFR within 11 Mpc is enhanced with respect to 
the cosmic SFR (shaded band) due to cosmic variance. Beyond 20 Mpc, the 
ineasurements and their mean, 0.017l^Q g[5g M0 Mpc^'* yr^^, are compat- 
ible with the cosinic SFR, and show no significant cosmic variance in this 
distance range. Shown for a Salpeter A IMF. 

agreement within uncertainties with the better-measured cos- 
mic SFR, except for the volume within < 10 Mpc which has 
a factor of < 2 enhancement. These results provide the im- 
portant basis on which we interpret the local SNR in Section 



3.3. Is the Fraction of Stars resulting in CC SNe too High? 

The predicted SNR depends on the number of stars that be- 
come CC SNe per unit stellar mass formed. To a first approx- 
imation, this depends on the mass range of stars that make 
CC SNe. There is observational evidence that a wide range 
of massive stars yield optically luminous CC SNe. Concep- 
tually, the mass range is controlled by two parameters, the 
lower mass limit that is the boundary between the formation 
of a white dwarf and a neutron star, and the upper mass limit 
that is less well defined but could be the boundary between 
forming a neutron star and a stellar mass black hole. This is 
only approximate, since stellar evolution and SN simulations 
indicate that parameters such as rotation and metallicity also 
affect the outcome of core collapse, and prompt black hole 
formation may be accompanied by some transient phenom- 
ena, even if dim. For example, a rare group of highly rotating 
stars that collapse to black holes f orm the current paradigm for 
gamma-ray burst central engines (IWoosleyll993h . We assume 
that parameters other than mass are averaged out in integrated 
populations and that the mass is the primary parameter for the 
CC SNe outcome. The existence of CC SNe progenitors of 
a given mass does not require or prove that all progenitors of 
that mass produce similar CC SNe. 

The lower mass limit is newly supported by direct obser- 
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vations of CC SNe progenitors. Using archival images, it 
has been possible to identi fy the stellar progenitors of some 
nearby CC SNe (see, e.g.. ISmarttll2009 [ for a review). Us- 
ing the stars' luminosities and other informati on, the stel- 
lar ma sses can be measured or limited. In ISmartt et al.l 
(l2009l) . the authors reviewed 20 progenitors of Type IIP SNe 
and statistically found the progenitor mass limit A/^^P — 
5M0. This is consistent with the highe st masses es- 
timated for white dwarf pr ogenitors, ~ 7 dKalirai et alj 
I2OO8; Willi ams et al.ll2009l) . Thus two different approaches 
seem to be converging to Afmin « 8 ± 1 Mq. This uncer- 
tainty affects the predicted SNR at the ±20% level. 

On the other hand, the upper mass limit is less certain. 
Fortunately, as long as the upper limit is large, it does not 
strongly affect the predicted SNR. For example, to fully ex- 
plain the normalization discrepancy in Figure [T] would require 
A'^max ~ 13 Mq for Mmin = 8 M©. Theoretically, the upper 
mass limit is Mmax ~ 40 M© for a sub-solar metallicity star 
dFrverll 19991) . Higher metalhcity generates stronger mass loss 
prior to collapse and predicts the formation of a Wo lf-Rayet 
(WR) star, which can extend A/max to 100 M© ( Hegeret al.1 
l2003h . The upper mass limit derived from the analysis of 



20 Type IIP SN p rogenitors is Af/JP^ = 16.5 ± 1.5 Mq 
dSmartt et al.ll2009h . No progenitors of Type Ibc SNe have 
been directly identified, but the spatial distribution of CC SNe 
shows that Type Ibc pr ogenitors must be more massive than 
those of Type IIP SNe dAnderson & Jamesll2008ll2009l) . It is 
widely expected that Type Ibc originate from evolved mas- 
sive WR stars that have shed their envelopes. From Galactic 
stellar clusters and stellar clusters in the LMC, such WR stars 
are estimated to hav e evolved f rom main-sequence masses 25 
Mq and above ( Massev et al.l ,2000, 200D- Very luminous 
Type Iln have been observed to arise from the core-collapse of 
very massive stars of the luminous blue variable (LB V) type 
jOal- Yam et al. 2007; Gal- Yam & Leonard 2009; SmithetaD 
I2OI lb . whose masses are > 20Mq. Recent studies of the pe- 
culiar SN 1961V support its nature as a true CC SNe, imply- 
ing its progenitor, a massi ve > 80 Mq star, produced an op- 
tically luminous CC SNe dKochanek et alj|2010l; ISmith et alj 
12010) . Therefore, there is substantial evidence for the upper 
mass limit being high. 

However, the outcome s of the mass range from 17 to 25- 
30 Mq remain uncertain. ISmartt et al.l (|2009.) noted that stars 
in this mass range are not found as CC SNe progenitors, and 
termed this the "red supergiant problem." Future progenitor 
studies may discover CC SNe progenitors in this mass range. 
For example, the progenitor of the recent Type IIL SN 200 9kr 
has been suggested to be ~ 20Mq dElias-Rosa et alj2010 l al- 
though Fraser et al. 2010a! find ~ 15Mq). It is also possible 
that the majority of these stars form bl ack holes at core col- 
lapse and lead to dim or dark CC SNe ( Smartt et a"Lll2009l) . In 
Section|5]we discuss how monitoring for the disappearance of 
high mass stars will allow obser vation of optically dark core 
collapses ( Kochanek et al.ll2008l) . 

In summary, we adopt the nominal mass range for opti- 
cally luminous CC SNe of 8-40 Mq, based on stellar and su- 
pernova simulations; the maximum mass range 8-100 Mq; 
and the conservative minimum mass range the combination 



of 8.5-16.5 Mq and 25-40 Mq. Compared to the nominal 8- 
40 Mq, the maximum and minimum mass ranges affect the 
predicted cosmic SNR by +10% and —30%, which are insuf- 
ficient to explain the supernova rate problem. 

3.4. Are Measurements Missing Luminous CC SNe? 

Luminous CC SNe (we henceforth define luminous CC 
SNe as CC SNe more luminous than —15 mag) that fall within 
the sensitivity of SN surveys can still be missed. Here, we dis- 
cuss the potential impact of these missed luminous CC SNe. 

Flux-limited surveys targeting a pre-selected sample of 
galaxies are naturally biased against small galaxies, and in- 
evitably result in some CC SNe being missed. Even the 
LOSS survey views only 45% of the total stellar ligh t ac- 
cessible from the survey position (Leaman et al.ll201ll) . In 
the early works of Cappellaro et al. (1999) and up dates by 
ICappellaro et al.1 d2005h and iBotticella et al.l (|2008|) . this has 
been corrected for by assuming that the CC SN rate per galaxy 
size (mass or luminosity) remains constant in galaxies having 
the same Hubble type or: B — K color Therefore, the volu- 
metric SNR could be determined even if the faint end of the 
galaxy luminosity function is not fully sampled, by measur- 
ing the CC SN rate in large galaxies and multiplying by an 
independently measured galaxy luminosity density that has 
been integrated to include small galaxies. However, more 
careful sampling by LOSS helped reveal the rate-size rela- 
tion, where the SN frequency is not linearly proportional to 
the galaxy size but rather size to the power of 0.4-0.6, the 
exact number depending on the SN Type and galaxy Hubble 
type dLi et al.ll201 lb). The rate-size relation is comparable 
to the trend of specific-SFR w ith galaxy mass dSalim et al] 
I2OO7I; ISchiminovich et al.ir2007i) and implies more emphasis 
on small galaxies. 

However, the rate-size relation has been measured with 
galaxies down to Lk « 2 x 10^° Lq {Mk ~ —22 A mag, 
approximately 1.3 mag fainter than the Milky Way), and ad- 
ditional effects in smaller galaxies may not addressed. In 
particular, there is mounting evidence that CC SNe in small 
galaxies are unlike those in large normal galaxies. Some 
types of rare CC SN e occur o nly in low-metallicit y dwarf 
galaxies dLe Floc'h et al. 2003; iFruchter et al.ll2006l) and re- 
sults from the Palomar Transient Factory (PTF), which rem- 
edy the bias against small galaxies by discovering the SNe 
first and the galaxy second, suggest systematic differences in 
CC SNe in small and normal galaxies. For example, more CC 
SNe have be en discovered in dwarf galaxies than predicted in 
lYoungetaL (2008) fro m the S FR and a CC SN mass range of 
8-50 Mq dArcavietai]|2010l). Also the d istributions of CC 
SNe types are different (l Arcavi et"ai1l2010l) . 

It is instructive to compare the distributions of CC SNe and 
SFR as functions of host galaxy properties. The fraction of 
CC SNe in dwarfs (defined as galaxies w ith Mr > — 18 ma g 
in the PTF) is 15/69 w 0.22 (Figure 3 of lArcavi et aPbOlOl) . 
even though the share of the total SFR in such galaxies is only 
0.1 (Figure 15 of Brinchmann et al. 2004). For a cut of AIr > 
— 19 mag, the fractions are 0.33 and 0.20, respectively. K 
corrections are insufficient to account for this difference, and 
the difference could be interpreted as an increased efficiency 
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of forming CC SNe in small galaxies. The increase is a factor 
of ^ 2 over that in normal galaxies, and the extra CC SNe 
would increase current SNR measurements. Since 20% of the 
SFR (in small galaxies) yields twice as many CC SNe, the 
increase in the SNR measurements is approximately 20%. As 
an extreme, we adopt an increase of 2 x 0.2/0.8 « 50%, i.e., 
assuming that all CC SNe in small galaxies are unaccounted 
for at present. A detailed comparison of the distributions of 
CC SNe and SFR requires more data and is beyond the scope 
of this paper The importance of small galaxies will become 
clearer with time, as SN surveys continue to discover more 
CC SNe in small galaxies. In particular, if the efficiency is 
enhanced over a wider range of galaxies, the increase in the 
SNR could be larger than discussed here. 

Finally, in all SN surveys, a short cadence is desirable for 
discovering SN during the most luminous epochs, as well as 
for obtaining a well sampled light curve for efficiency calcu- 
lation purposes. The cadences of surveys, shown in Table [T] 
are generally shorter than the peak durations of CC SNe. For 
example. Type Ibc typically drop by a magnitude in less than 
1 month, while Type IIP evolve more slowly, with a plateau 
remaining within 1 mag of the peak for about 2-3 months. 

In summary, incomplete surveying of small galaxies may 
lead to SNR measurements underestimating the true SNR. 
Based on current data on dwarf galaxies, the effect may be 
a 20-50% increase in the SNR measurements. The impor- 
tance of small galaxies remains to be investigated further with 
more data, including its caveats (Section l4!4] i. 

3.5. Are CC SNe Dust Corrections Insufficient? 

Host dust extinction makes CC SNe appear dimmer and 
correcting for missing obscured CC SNe is often the most 
uncertain ingredient in SNR measurements. Here we review 
dust corrections performed in SNR measurements. 

The dust correction applied to SNR measurements ranges 
from a few tens of percent at low redshifts to a factor of 2 
at high redshifts. A ccording to the semi-analytic model of 
iHatano et al] (Il998l) . the distribution of CC SN extinction due 
to host galaxy inclination peaks at low extinction but has an 
extended tail reaching many magnitudes (19.8 mag for a per- 
fectly edge-on galaxy). This model was applied in SNLS, 
who found that it in creased their m e asured SNR by « 15% 
dBazin et alj|2009h . iBotticella et al. ' ('2008') used a statistical 
approach similar to|Hatano et al. ( 1998) and found an increase 
of a factor of ^ 2. Although the authors also present an ex- 
treme dust scheme which would in fact find agreement with 
the SNR predicted from the cosmic SFR, they point out that 
such an extreme dust scenario is not strongly mot ivated and is 
meant as an upper limit only. At higher redshif t . iDahlen et all 
(|2004 combined the model of iHatano et all ( Il998h with a 
starburst extinction law, reporting increases of 60% and a fac- 
tor of ^ 2 for their two measurements at z ~ 0.3 and z « 0.7, 
respectively. 

As pointed out in 'Ma nnucci et alj (|2007'), at high redshift 
there is additional extinction due to starburst galaxies and 
highly star-forming galaxies (luminous and Ultraluminous 
IR galaxies, which we collectively call ULIRG). Identified 
by their strong FIR emission and high SFR, ULIRG should 



house many more CC SNe per galaxy, but we only detect a 
small fraction of CC SNe because of the higher dust obscura- 
tion. It has been estimated that as many as 60%-90% of the 
CC SNe could go undiscovered (Mannucci et al. 200^, al- 
though m ulti-wavelength studi es are starting to discover some 
of these jKankare et al.l 12008*). Because the importance of 
ULIRG increases with redshift, the fraction of CC SNe that 
are expected to be missing increases from ~ 5% at z « to 
20%-40% at z « 1 ( Mannucci etallbOOTl) . 

In summary, correction for host galaxy obscuration remains 
uncertain. Extreme corrections that would explain the super- 
nova rate problem are not prohibited, in particular at high red- 
shifts where the excess extinction due to starburst galaxies is 
considered. However, this does not easily explain the super- 
nova rate problem at low redshifts, where starburst galaxies 
provide only a small fraction of the SFR. We discuss missing 
dim CC SNe in more detail next. 

3.6. Is the Contribution from Dim CC SNe Significant? 

Type IIP SNe are the most common type of CC SNe 
and also the most varied in terms of luminosity (e.g., 
iRichardson et al.ll2002t lHamuvll2003l) . The luminosity of a 
typical Type IIP SN remains nearly constant for a relatively 
long duration of ~ 100 days (the plateau phase), after which 
it dro ps sharply, mark ing the transition to the nebular phase 
(e.g., Patat et al]|l994 ^. It is expected that the radius of the 
progenitor plays a key role in shaping both the plateau length 
and luminosity (lArnettlll980l) . For example, the unusually 
dim ligh t curve of SN 1987 A is attributed to its compact pro- 
genitor (lArnettl lT987). Recently, SN searches have revealed 
a class of Type IIP SNe that have peak and t ail luminosities 
that are even lower than tho se of SN 1987 A (IZampieri et al.l 
I2OO3I: IPastorello et"aDl2004 . The prototype of these ^^^Ni- 
poor (< 0.01 M©) lo w-velocity Type IIP SNe is SN 1997D 
(My w -14.7 mag: lTuratto et al]ll998h . an d there are now 
more t han a d ozen similar CC SN e (see, e.g., IPastorello et al.l 
'2004' '2006t iFraser et al.l l2010bl). The dimm est are SN 
1999br (Mr w -13.5 mag: IPastorello et al.l 12004) and 
PTFlOvdl/SN 2010id (Mb. « -14.0 mag: iGal-Yam et alJ 
l20Tll) . although even dimmer suspects exist, e.g., the tran- 
sient in M85 (Mb. ^ —12 mag: Pastorello et al. 2007, see also 
iKulkarni et alj20()7l for an alternative scenario). Additionally, 
CC SNe dimmer than Mb ~ —13.5 mag have been inferred 
from the non-observation of SNe associated with several 
low-r edshift long-GRBs ( Fvnbo et al.H2006l: [Pella Valle et alj 
I2OO6I) . There is a population of dim (My < — 12 mag) Type 
Iln SNe, but these are su spected to be LBV outbursts and not 
core collapses (e.g., Maund et al.ll2006l:ISmith et al.ll20I0l) . 

In addition to the cadence and galaxy completeness effects 
mentioned in Section|33| SN surveys are biased against low- 
luminosity CC SNe like SN 1997D. The Hmiting absolute 
magnitudes of SN surveys are typically M w —15 to —16 
(Table [1]). Although SN surveys can discover dimmer CC 
SNe, the rapidly falling detection efficiency and the smaller 
volume from which dim CC SNe (we define CC SNe that are 
dimmer than M = — 15 mag as "dim" CC SNe) may be confi- 
dently discovered makes collecting dim CC SNe challenging. 
To correct for missing dim CC SNe, a SN luminosity func- 
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Fig. 3. — Differential catalog SNR density as a function of distance. The 
total, luminous, and dim catalog SNR are shown by points (lowest distance 
bin has been shifted slightly in distance for clarity). The total catalog SNR 
is larger than the sum of the luminous and dim catalog SNR (see the text). 
Vertical error bars are statistical only. The SNR predicted from the cosmic 
SFR and directly measured are shown and labeled. The total catalog SNR is 
comparable to the upper band out to ~ 25 Mpc. The luminous catalog SNR 
is comparable to the lower band out to ~ 25 Mpc. The dim catalog SNR 
is large out to ~ 10 Mpc, indicating a significant fraction of dim CC SNe, 
independent of the absolute SFR enhancement within ~ 10 Mpc. 

tion that is complete to dimmer CC SNe must be adopted. 
Here, we revisit the importance of dim CC SNe. We make 
use of SN discoveries recorded in SN catalogs. The SN cata- 
log combines results from SN surveys and amateur discover- 
ies, resulting in a non-uniform sample of CC SNe of various 
qualities, bands and uncertainties. However, it is illustrative 
because it sets a strong lower bound on dim CC SNe. We first 
discuss the local (< 100 Mpc) volume for a qualitative anal- 
ysis in Section [3.6.1l Then, we focus on the very local (< 10 
Mpc) volume for a more quantitative analysis in Section l3.6.2l 

3.6.1. How Ubiquitous are Dim CC SNe? 

We start wi th the Sternberg Astr onomical Institute (SAl) 
SN Catalog dBartunov et al.n2007b . selecting CC SNe of 
Types IIP, IIL, lln, lib, and Ibc. We select only the most 
recent decade, from 2000 to 2009, sin ce previous decades 
have severe incompleteness problems dHoriuchi & BeacomI 
|201(^. We cross check the classification with the Harvard SN 
databas^IB. which results in tagging some unclassified SNe as 
CC SNe. We update magnitudes to peak magnitudes where 
possible by investigating the literature and circulars^*!. We 
update SN host galaxy distances with redshift-independent 
distanc e measures from the Extragalactic Distance Database 
(EDD: ITuUv et alj|2009l) . When these are not available, we 
used distance estimates from NEeQ The uncertainties in the 
distances, which can reach ±10%-20%, affect the absolute 

http://www.cfa.harvard.edu/iau/lists/Supemovae.html 
' ' http://www.cfa.harvard.edu/iau/cbat.html 
' ^ http ://nedw w w. ipac . caltech.edu/ 



magnitudes. As we discuss in Section 14.31 this does not 
strongly affect our results. 

From the SAl catalog we discard CC SNe that seem to 
clearly be LBV transients (e.g., SN 2000ch, SN 2001ac, SN 
2002kg, SN 2003gm, SN 2006fp, SN 2007sv, SN 2009ip). 
On the other hand, we include the possible SN impostors 
SN 2008S, NGC300-OT, and SN 2002bu, whose nature as ei- 
ther true CC SNe or extreme stellar outbursts remains subject 
to debate; we discuss the alternative assumption in Section 
14.41 Next, we check the discovery phases of CC SNe dimmer 
than —16 mag to identify CC SNe that are dim due to being 
discovered post-peak. For this identification we adopt the fol- 
lowing criteria: SN IIP discovered more than 2 months post- 
peak, SN IIL over 1.5 months, and all SN lln, lib, Ibc over 1 
month. These are the typical time scales over which t he nom- 
inal lig ht curves remain within 1 mag of the peak (iLi et al.l 
12011 at) . We identify 36 CC SNe that fit this category during 
2000-2009 within 100 Mpc. 

From the modified SAl sample we construct three sets of 
"catalog" SNR. The first is the definitely luminous catalog 
SNR. For this we select CC SNe with absolute magnitudes 
below —15 mag. The second set is the definitely dim catalog 
SNR. For this we select CC SNe above —15 mag. We do not 
include dim CC SNe that have been discovered in their late 
phases, as the true peak magnitude is uncertain. The third 
and final set is the total catalog SNR. This includes all CC 
SNe, including definitely luminous CC SNe, definitely dim 
CC SNe, and late-time discoveries. Thus, the total catalog 
SNR is larger than the sum of the definitely luminous and 
definitely dim catalog SNR. 

All catalog SNR should be considered lower limits since 
they are derived from simple counting of likely incomplete 
SN discove ries. In addition to increas ing incompleteness in 
older data (iHoriuchi & BeacoinlbOlOl) . there is a strong ob- 
servational bias against CC SNe in the southern hemisphere. 
In the latest decade and within 100 Mpc, there are almost 1.9 
times more CC SNe discovered in the Northern hemisphere 
than the Southern hemisphere. There is also a strong bias 
against CC SNe in small galaxies. Considerations of the ob- 
servational efficiency and coverage will only increase the cat- 
alog SNR. 

If the SN discoveries are sufficiently complete, the catalog 
SNR should be flat with distance. This is because they are 
number densities, and because the local SFR did not reveal 
significant cosmic variance, except in the very local volume 
within ^ 10 Mpc (Section [3.2l ). In other words, any signifi- 
cant decrease of the catalog SNR at any distance signifies an 
incompleteness of the SN discoveries at that distance. 

In Figure |3] the catalog SNR are binned in distance and 
compared to the cosmic SNR predicted from the SFR and di- 
rectly measured. Neither host nor Galactic extinction correc- 
tions have been applied to the catalog SNR. Correcting would 
shift a small number of dim CC SNe to luminous CC SNe, but 
would not change qualitative conclusions; in Section D.6.21 
where we discuss quantitative results, we apply Galactic ex- 
tinction corrections. Each bin contains 20 luminous CC SNe, 
except for the smallest-distance bin, which contains 10 lumi- 
nous CC SNe. 
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TABLE 2 

Very Local CC SNe During 2000 to 2009 Inclusive 



SN 


Galaxy 


Type 


D (Mpc) 


E{B - V) 


Absolute Magnitude" 


Discovery Phase 


References 




SN 2002bu 


NGC4242 


Iln" 


5.8 


0.012 


Mr ^ 


-14.1 


Early 


Homoch (2002); Smith et al. (2010) 


SN 2002hh 


NGC6946 


IIP 


5.9 


0.342 


Mr ^ 


-14.3 


Early 


Pozzo et al. (2006) 


SN 2002kg 


NGC2403 


LBV 


3.2 


0.04 


Mv ~ 


-9 


Not CC SN 


Maund et al. (2006) 


SN 2004am 


NGC3034 (M82) 


IIP 


3.5 


0.159 


Mr < 


-12.2 


3 months 


Singer et al. (2004) 


SN 2004dj 


NGC2403 


IIP 


3.2 


0.04 


Mr ~ 


-16.0 


1 month 


Vinko et al. (2006): Zhans et al. (2006) 


SN 2004et 


NGC6946 


np 


5.9 


0.342 


Mr ^ 


-17.6 


Early 


Sahu et al. (2006) 


SN 2005af 


NGC4945 


IIP 


3.6 


0.177 


Mr ~ 


-15.4 


1 month 


Jacques & Pimentel (2005) 


SN 2005at 


NGC6744 


Ic 


7.1 


0.043 


Mr ~ 


-15.1 


2 weeks 


Martin et al. (2005b) 


SN 2008bk 


NGC7793 


IIP 


4.1 


0.019 


Mr- 


-15.5 


1 month 


Monard (2008a); Morrell & Stiitzineer (2008) 


SN 2008iz 


NGC3034 (M82) 


II? 


3.5 


0.159 


no optical 


Radio only 


Bmnthaler et al. (2010) 


SN 2008S 


NGC6946 


Iln'' 


5.9 


0.342 


Mr S:! 


-13.3 


Early 


Botticella et al. (2009); Smith et al. (2009) 


NGC300-OT 


NGC300 


Iln'' 


1.9 


0.013 


Mv ~ 


12.3 


1 month 


Bond et al. (200^) 


SN 2002ap 


NGC0628 


IcPec 


9.0 


0.07 


Mr ^ 


-17.8 


Early 


Gal- Yam et al. (2002) 


SN 2003gd 


NGC0628 


IIP 


9.0 


0.07 


Mr ~ 


-16.7 


2 months 


McNaught (2003) 


SN 2005cs 


NGC5194(M51) 


IIP 


8.4 


0.035 


Mr 


-15.4 


1 month 


Pastorello et al. (2006) 


SN 2007gr 


NGC1058 


Ic 


9.9 


0.062 


Mr ^ 


-17.4 


Early 


Valenti et al. (2008) 


SN 2008ax 


NGC4490 


lib 


9.6 


0.022 


Mr ^ 


-16.6 


2 weeks 


Pastorello et al. (2008) 


SN 2009hd 


NGC3627 (M66) 


IIP 


8.3 


0.032 


Mr ^ 


13.9 


Early 


Monard (2009) 


SN 2001ig 


NGC7424 


lib 


11.5 


0.011 


Mr ^ 


-17.3 


Early 


Bembiick et al. (2002) 


SN 2003ie 


NGC4051 


II 


12.2 


0.013 


Mr < 


-15.6 


Uncertain 


Arbour & Boles (200J) 


SN 2003jg 


NGC2997 


Ibc 


11.3 


0.109 


Mr ~ 


-14.1 


Few weeks 


Bigss et al. (200^ 


SN 2007it 


NGC5530 


IIP 


11.7 


0.116 


Mv ~ 


-18.7 


early 


Poimanski (2007) 


SN 2008eh 


NGC2997 


Ibc? 


11.3 


0.109 


Mr ~ 


-15.3 


1 month 


Monard (2008b) 


SN 2009ib 


NGC1559 


IIP 


12.6 


0.03 


Mr ^ 


15.9 


Early 


Pignata (2008) 



Note. — CC SNe that are definitely dimmer than —16 mag are shown in bold type. CC SNe with magnitude limits may be dim or luminous CC SNe, but cannot 
be convincingly identified with the data available. The list has been divided into blocks by distance (7.9, 10, and 12.6 Mpc). The first block and second block contain 
equal volumes, while the third block contains twice that. 

" Shows the peak magnitudes for pre-peak discoveries, and the discovery or inferred peak magnitude for others. Values have been corrected for Galactic extinction. 
Unfiltered images are categorized as i?,-band. 
Possible SN impostors, which may be CC SNe or stellar outbursts. 



The luminous catalog SNR should be compared to the mea- 
sured cosmic SNR. The luminous catalog SNR is generally 
decreasing with distance but it is flat out to ^ 25 Mpc, sug- 
gesting reasonable completeness to that distance. The nor- 
malization of the flat section is comparable to the measured 
cosmic SNR on larger scales by LOSS, supporting this inter- 
pretation. There are 50 luminous CC SNe in the first 3 dis- 
tance bins. The smallest-distance bin shows an enhancement 
as expected from the SFR. 

The total catalog SNR should be compared to the predicted 
cosmic SNR. The total catalog SNR is flat out to ~ 25 Mpc 
and shows a normalization that is only slightly lower than 
expectations from the cosmic SFR. This supports our ear- 
lier c laims that the true cosmic SNR is a s large as expected 
(e.g.. lHoriuchi et al.ll2009t lBeacomll2010l) . There are 75 CC 
SNe within the first three distance bins. Again, the smallest- 
distance bin shows an enhancement with respect to the pre- 
dicted cosmic SNR, as expected from the SFR. 

Finally, the dim catalog SNR falls with distance at all dis- 
tances, a sign of incompleteness. At the smallest-distance bin, 
where dim CC SNe are least likely to be missed, the dim cat- 
alog SNR is, surprisingly, just as large as the luminous SNR 
measured by LOSS. There are 4 dim CC SNe in this bin. This 
normalization is large enough to help solve the supernova rate 
problem, at least at z « 0, although the enhancement of the 
SFR in this distance range needs to be taken into account. 

Catalog CC SNe and Type la supernovae (SNe la) can be 



used to make an inde pendent test of the no rmalization of the 
cosmic CC SNe rate jHoriuchi et al.ll2009l) . SNe la are lumi- 
nous and have been intensively sought. Their cosmic rate his- 
tory is better measured than that of CC SNe, and the distance 
out to which their catalog is apparently complete is larger 
dHoriuchi & BeacomI 12010) . We consider the SN la to CC 
SN ratio, comparing rates from cosmic surveys near z « 
(for larger z, the effect of SN la progenitor delays changes the 
ratio) and SN catalogs. From our fits to the measured cosmic 
rates, we find a SN la to CC SNe ratio at z « of 0.18. 
For comparison, the same ratio for the LOSS results is ~ 0.43. 
This difference of a factor of ^ 2 in the ratios is another state- 
ment of the supernova rate problem. The same ratio for the 
SN catalog should largely cancel incompleteness or local en- 
hancement effects and is independent of the assumed cosmic 
SFR. We use a maximum distance large enough to contain 
enough SNe la and small enough that the CC SNe are rea- 
sonably complete. For SNe within 20 Mpc, we find a ratio 
of 0.25; this increases to 0.30 by 25 Mpc and rises to 0.35 
at larger distances. We interpret this as confirming that cos- 
mic CC SNe are as common as predicted and that many dim 
ones are missed at increasing distances in catalogs and cosmic 
surveys. 

In summary, we analyze the local (< 100 Mpc) volume 
to study dim CC SNe. The SNR follows the trend expected 
from the SFR out to 25 Mpc. We find that dim CC SNe are 
severely incomplete. Despite this, dim CC SNe are numerous 
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TABLE 3 

Summary OF Very Local CC SNe 4 



Number of CC SNe 0-7.9 Mpc 7.9-10 Mpc 10-12.6 Mpc 



Expected total (cosmic rate) 2.9 ± 0.5 2.9 ±0.5 5.8 ± 1 

Expected (LOSS rate) 1.5 ±0.4 1.5 ± 0.4 3.0 ± 0.9 

Total observed CC SNe 12" 6 6 

M < -16 2 4 2 

-16 < Af < -15 3 1 2 

-15 < M (dim) 2 1 1 

Possible SN impostor 3 

Poor-quality CC SNe 2 1 



p-i— I— I— p-i— I— I— p-i— I— I— p-i— I— I— p-i— I— i-q 




Dim fraction f^. 

dim 

Fig. 4. — Probability density distribution for the dim fraction (/dim), given 
the observed number of total and dim CC SNe, for three volume cuts. Shown 
for CC SNe with poor-quality CC SNe conservatively removed; see the text 
and Table[3] The dim fraction is significantly higher than previous estimates 
from larger distances (vertical bands). Top: for CC SNe dimmer than —16 
mag. Bottom: for those dimmer than —15 mag. 

light curve most closely resembles a dim Type Ibc CC SNe 
rather than a SN 2002cx-like peculiar Type la. 

It is not feasible to estimate the peak magnitudes of some 
CC SNe. For example, SN 2004am was discovered in the neb- 
ular phase, and the discovery phase of SN 2003ie is unknown. 
We collectively treat these CC SNe as poor-quality CC SNe. 
We exclude these CC SNe from estimates of /dim, since we 
want to estimate the dim fraction from a clean sample of SNe. 

Within 10 Mpc, there are 18 CC SNe including 1 LBV 
(SN 2002kg). Excluding poor-quality CC SNe leaves 16 
CC SNe, of which 6 are more luminous than —16 mag 
(SN 2002ap, SN 2003gd, SN 2004dj, SN 2004et, SN 2007gr, 
SN 2008ax) and 10 are dimmer than —16 mag. Of the latter, 

3 are highly reddened (SN 2002hh, SN 2008iz, SN 2009hd), 

4 are intrinsically dim (SN 2005af, SN 2005at, SN 2005cs, 
SN 2008bk), and 3 are possible SN impostors (SN 2002bu, 
SN 2008S, NGC300-OT). The number of CC SNe dimmer 
than —15 mag is 6, so the dim fraction is />-i5 = 6/16 
(/>-i6 — 10/16). If we include poor-quality CC SNe and 
allow them to be either luminous or dim CC SNe, the dim 
fraction varies between />_i5 =6/18 and />_i5 = 8/18. 

We expect the dim fraction to decrease with distance, since 
dim CC SNe are more easily missed at larger distances. Go- 
ing out to twice the volume (i.e., within 12.6 Mpc) we find 6 
more CC SNe (5 after excluding poor-quality CC SNe). There 
is 1 dim CC SNe (SN 2003jg), bringing the dim fraction to 
/>-i5 = 7/21 (/>_i6 = 13/21). Considering half the vol- 
ume (i.e., CC SNe within 7.9 Mpc) we find 12 CC SNe (10 
after excluding poor-quality CC SNe), with />-i5 = 5/10 
(/>-i6 ~ 8/10). The CC SNe tallies are summarized in Ta- 
ble [3 



Note. — Expected number of CC SNe are shown for comparison. They 
are derived from measurements at larger distances, and should be multiplied 
by < 2 within ~ 10 Mpc to reflect the overdensity of SFR at these distances. 
The uncertainty is from the 1 a uncertainty on the cosmic SFR or the uncer- 
tainty on the SNR measurements. 
" 9 CC SNe if possible SN impostors are not included. 

in the very local (< 10 Mpc) volume. In the next Section, we 
discuss the fraction of dim CC SNe in more detail. 

3.6.2. What is the Fraction of Dim CC SNe within ~ 10 Mpc? 

We use the very local (< 10 Mpc) volume to quantify the 
importance of dim CC SNe. It is not feasible to extend to 
greater distances because of the severe incompleteness of dim 
CC SNe. We choose a nominal distance of 10 Mpc, which 
yields sufficient statistics, and discuss how halving and dou- 
bling the volume changes the results. As discussed in the 
previous section, the normalization of the dim catalog SNR 
within 10 Mpc is high. However, the very local volume has an 
enhanced absolute SFR. Therefore, we conservatively use the 
fraction of dim CC SNe, /dim = Ndim/{Ndim + A^iuminous), 
the ratio of definitely dim CC SNe over the sum of definitely 
dim and definitely luminous CC SNe. Here, the subscript dim 
refers to objects fainter than dim-th magnitude. We will focus 
on a dim- luminous magnitude cut of —15 mag, but also show 
— 16 mag for illustration. 

The very local CC SNe are shown in Table |2] Here, the 
magnitudes ar e corre cted for Galactic extinction fol l owing 
ISchlegel et all (Il998h and assuming a I Cardelli et al.1 (Il989h 
Galactic reddening law. Distance estimates for CC SNe host 
galaxies vary among sources, and we adopt the estimates used 
in dedicated papers on the CC SNe when this is available. 
Otherwise, we continue to use the distances from FDD or 
NED as explained in Section [J!6l 

Several CC SNe are dim due to being highly extinguished 
by host dust, e.g., SN 2002h h with Ay » 5 mag and SN 
2008ax with ~ 1 mag ( Pozzo et alJl2006l: iMattila et alJ 
12004) . The extreme case is SN 2008iz, which was discovered 
in radio observations and is not seen in the optical, suggesting 
an explosion behind a large gas o r dust cloud in the central 
part of the host galaxy (Brunthaler et al.ll2010l) . We assume 
SN 2008eh to have been a CC SN even though it lacks a spec- 
troscopic classification. Despite its relatively early discovery 
(it was not observed in the field 1 month prior to discovery), 
it was dim (M^ ^ -15.3: iMonar d 2008b). From archival 
Spitzer images of its host galaxy NGC 2997, we locate the 
SN position in a spiral arm, next to a star forming region. Its 



12 



Horiuchi et al. 



In Figure|4] we plot the probability density function of /dim 
given the observed number of total and dim CC SNe, for dif- 
ferent distance cuts. We use the binomial distribution in A^dim 
for a fixed iVtot and varying /dim, using a flat prior in /dim 
over the interval 0-1. On the top panel, we show a magni- 
tude cut of —16 mag for illustration. Due to the broad distri- 
bution of /dim, we are not able to derive a precise value for 
/dim with the present data: the fractions are />-i5 ~ 30%- 
50% (/>-i6 60%-80%), but lower and higher values are 
not improbable. However, our estimates are clearly larger 
than previous studies based on more systematic samples at 
higher distances which were not as sensitive to the dim end of 
the CC SNe luminosity function (/>-i5 3%-20%; Sec- 
tion |43]l. Our estimates decrease to />-i5 ^ 20%-40% 
(/>-i6 ^ 50%-70%) if possible SN impostors are excluded 
from our list of CC SN. This is still larger than previous stud- 
ies. 

In summary, the very local 10 Mpc volume shows the dim 
fraction could be as high as />-i5 ^ 50%. This is signif- 
icantly higher than previous estimates, which range between 
/>-i5 ^ 3%-20%. Thus the SNR measurements are scaled 
up by (1 - 0.2)/(l - 0.5) = 1.6 (a factor of 1.9 if the orig- 
inal dim fraction is 5%), largely solving the supernova rate 
problem. While the very local 10 Mpc volume is special in 
its absolute SFR, this should not affect the fraction of dim CC 
SNe. 

4. DISCUSSION 

In Figure |5] we show the supernova rate problem at z w 0. 
The bands reflect the nominal uncertainties on each quantity: 
for the predicted SNR it is the la uncertainty in the cosmic 
SFR fit by Hopkins & Beacom ( 2006), while for the measured 
SNR it is the combined statistical and systematic uncertainty 
of the fit to SNR measurements. Although the supernova rate 
problem remains over the entire redshift range where SNR 
measurements are available, we focus here on the z « range 
where there is most data. 

4.1. Summary of Potential Explanations 

First, to investigate whether the predicted SNR is too large, 
we explore the intrinsic uncertainties in the SFR measure- 
ments (IMF, obscuration, and conversion factor) and the addi- 
tional factors when the SFR is converted to a SNR (IMF and 
CC SN mass range). The results are shown in Figure |5]by 
dashed arrows. 

• IMF. We vary the IMF shape from the Salpeter A 
IMF to a classic Salpeter IMF and a flatter Baldry- 
Galzebrook IMF. This only has a small effect on the 
predicted SNR. 

• Obscuration. Obscuration correction is redshift depen- 
dent. Our adopted SFR compilation is corrected for 
obscuration by adding the FIR-derived SFR to the UV- 
derived SFR. We consider uncertainty due to two poten- 
tial FIR contaminants: first, from AG N, which has bee n 
estimated to be < 10% of the FIR jSilva et alj|2004 . 
and second, from cirrus emission, which we assume 
to be < 50% of the FIR (.Lonsdale Persson & Heloi] 
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Fig. 5. — Measured and predicted SNR at 2: 0. Shaded bands reflect un- 
certainties: statistical and systematic combined for the measured SNR, and 
the Icr uncertainty in the SFR for the predicted SNR. The results of varying 
the inputs for the predicted SNR ai'e shown by dashed arrows, while the pos- 
sible effects of missing CC SNe are shown for the measured SNR by solid 
arrows. Stacked arrows reflect different assumptions (see the text). The con- 
tribution from dim CC SNe can potentially increase the measured SNR to 
bridge the normalization discrepancy, although other effects can also be im- 
portant. 

Il987h . We thus assume a reduction to the predicted 
SNR of 25%. Extinction uncertainties are unlikely to 
reduce the SNR prediction normalization by a full fac- 
tor of 2 (see Section [TTI for integral constraints on the 
low redshift SFR normaUzation). 

• SFR conversion factor. We assume formal uncertainties 
of ±15%. Although conversion uncertainties for young 
star forming galaxies are larger, such galaxies are rare 
in the low-z universe dSoifer & Neugebauerlll99l]) and 
can be safely neglected. 

• CC SN mass range. The nominal mass range is 8- 
40 Mq. We use a minimum mass range taken from 
observations of CC SN progenitors: 8.5-16.5 M0 and 
25-40 M0. We stress that this range does not account 
for the progenitors of all CC SNe Types. Further SN 
observations are needed to clarify this situation. We 
adopt the maximum mass range of 8-100 M©. 

Therefore, the uncertainties affecting the predicted SNR are 
generally smaller than the normalization discrepancy. 

Second, to investigate whether the measured SNR is too 
small, we explore the contribution from missing CC SNe (in- 
complete galaxy sample and dim CC SNe). This is shown by 
solid arrows in Figure |5] 

• Incomplete galaxy sample. This only affects SN 
surveys targeting a pre-selected sample of galaxies. 
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Though careful work is done to correct for the bias 
against CC SNe in small galaxies, recent results suggest 
CC SNe in dwarf galaxies are unlike those in normal 
galaxies. We adopt a tentative factor of two enhance- 
ment in the CC SNe formation efficiency in dwarfs, 
which would yield a 20%-50% increase in the total 
SNR, as shown by stacked arrows. Detailed studies 
need to await more data. 

• Dim CC SNe. There is an observational bias against 
the discovery and correction of dim CC SNe. In the 
very local (< 10 Mpc) volume, where dim CC SNe are 
least likely to be missed, the fraction of dim CC SNe 
is higher than in SN surveys targeting larger distances. 
Using the yeai's 2000-2009, we find />_i5 ~ 30-50% 
(if possible SN impostors are excluded from our list of 
CC SNe, we find 20%-40%), whereas SN surveys find 
3%-20%. Therefore, we scale up the SNR measure- 
ment accordingly by a factor of between (1 — 0.2) /(I — 
0.3) w 1.14 and (l-0.03)/(l- 0.5) « 1.9, shown by 
stacked arrows. 

Therefore, if possible SN impostors are true CC SNe, the frac- 
tion of dim CC SNe can be sufficiently high that dim CC SNe 
largely explain the supernova rate problem. If possible SN 
impostors are not true CC SNe, /dim is smaller, and dim CC 
SNe do not fully bridge the normalization discrepancy. 

In Figure|5]we show how the different effects compare. The 
supernova rate problem is likely explained by a combination 
of the effects shown. 

4.2. How Representative is the Very Local Volume? 

The validity of using the dim fraction derived from very 
local (< 10 Mpc) CC SNe relies on how representative the 
very local volume is. We discuss several issues. 

First, distance uncertainties for the very local galaxies can 
be up to ±25% in some cases. This has two effects: (1) the 
absolute magnitude is affected and can cause erroneous cate- 
gorization of CC SNe near the dim-luminous magnitude cut 
and (2) CC SNe can be categorized in the wrong distance bin. 
We have considered the dim fraction using different magni- 
tude cuts and different distance cuts, showing that it is consis- 
tently higher than previous estimates. 

Second, the very local volume shows an enhancement in its 
SNR (Figure[3]and Table|3j. The enhancement is significant: 
within 12.6 Mpc we observe 14 luminous CC SNe while the 
expectation is 6, an ^ 0.2% occurrence. The enhancement 
is a factor of ^ 2, which is comparable to the factor of 1.7 
enhancement observed in the SFR in similar distances. The 
enhancement is observed not only in the luminous CC SNe 
but also the total CC SNe: there are 23 discovered (including 
3 possible SN impostors) while the expectation is 12. This 
can be interpreted as dim CC SNe are not severely incomplete 
within these distances. 

Finally, to apply the dim fraction at cosmological dis- 
tance requires some care. The IMF at high redshift remains 
poorl y constrained (iNagashima et al.ll2005t iLe DelUou et al.l 
l2006l) . and studies find t hat IMF ev olution may be required 
at z > 0.5 (e.g.. .Wilkins et alj|20"08a) . If dim CC SNe are 



mass-dependent, an evolving IMF will lead to an evolving 
/dim- Also, redshift dependence may arise from metallicity 
dependence. A more important issue is obscuration by dust, 
which is known to be redshift-dependent (Section [3.5b . Our 
work in the very local volume is mainly on normal galax- 
ies, whose importance is overtaken by ULIRG at high red- 
shifts (Mannucci et al. 2007). The higher dust obscuration in 
ULIRG causes as man y as 90% of the CC S Ne to be missed 
by current SN surveys ( Mannucci et al.ll2003l) . Assuming that 
the dim fraction in ULIRG is 90%, the weighted dim fraction 
increases from ~ 50% at z « to ~ 70% at 2; ~ 1. 

In summary, we find that the very local (< 10 Mpc) vol- 
ume has a high absolute SNR, which is related to an enhanced 
SFR. The dim fraction should be independent of the absolute 
normalization, and thus our estimate may be used at low red- 
shifts. However, the dim fraction may increase with redshift. 

4.3. Dim CC SNe: Are They Compatible with Other Studies? 

One of the earliest discussions of dim CC SNe was based 
on historical CC SNe in the Milky Way and nearby galaxies 
(ISchaefeilll996l) . It was suggested that the CC SN luminosity 
function rises toward dimmer objects, with perhaps a third or 
more of the CC SNe dimmer than ]\Ib = —15. However, the 
completeness of historical SNe has been difficult to assess, 
and the SN sample contains significant extinction and type 
uncertainties. Some of the dim CC SNe are now known to 
have been LBV, e. g., SN 1954J, for wh ich the surviving star 
has been observed jVan Dvk et al ] l2005b . 

In more recent st udies, the fractions are smaller 
[Richardson et al.l (|2002|) report the mean and standard devia- 
tion of Type IIP, Type IIL, and Type Ibc SNe luminosity func- 
tions to be Mb « -16.18 (a « 1.23), Mb ~ -17.49 (cr w 
1.29), and Mb « -17.37 (a « 0.88) (adjusted to our chosen 
cosmology). Using the type fractions of the volumetric SN 
sample of LOSS, the fraction of C C SNe dimmer than M b — 
-15 (-16) is only 12% (34%). In lBotticella eTal] ilOOA the 
authors adopt somewhat narrower luminosity functions which 
result in a smaller dim fraction of 3% (13%). The SNLS uses 
an absolute ma gnitude scale rela tive to the mean magnitude of 
Type la SNe (iBazin et al.ll2009 l). Adopting My = -19 mag 
as the mean Type la SN peak magnitude, the fraction of CC 
SNe dimmer than My = -15 (-16) is 11% (30%). Here we 
use their luminosity function that has been corrected for the 
smaller target volume for dim CC SNe (data points of their 
Figure 12). The recently published luminosity functions of 
LOSS contain large numbers of dim CC SNe, including sev- 
eral in TableH e.g., SN 2002hh, SN 2004am, SN 2007gr, and 
SN 2008ax. From their luminosity functions, the fraction of 
CC SNe dimmerthan Mr = -15 (-16) is 20% (40%). Here, 
we have summed their luminosity histograms, rather than in- 
tegrate their Ga ussian fits which the authors caution are poor 
fits to the data ihi et al.ll201 Ibl) . In summary, the dim frac- 
tion in the hterature is />_i5 3%-20% (/>-i6 ^ 10%- 
40%). In Figure|4]we label the fractions according to SNLS 
and LOSS, the two most recent and reliable measurements. 

Using the very local SN catalog data, we estimate the dim 
fraction to be />_i5 - 30%-50% (/>_i6 - 60%-80%), 
larger than those of previous studies. Furthermore, we have 
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assessed the completeness of the catalog CC SNe (Section 
I3.6.11 I. The likelihood of our high /dim being a simple statis- 
tical fluctuation is small given the sample size, e.g., the like- 
lihood of fluctuating />_i5 - 3% (20%) to />_i5 > 50% in 
a sample size of 10 CC SNe is 0.0005% (3%). If possible SN 
impostors are excluded, f>-i5 ^ 20%-40% (/>-i6 ^ 50%- 
70%). This is in better agreement although still somewhat 
larger than the most recent luminosity function of LOSS. It 
is still significantly larger compared to other luminosity func- 
tions. 

An important reason is the small distance cut adopted. In 
SN surveys targeting larger distances, the dimmest CC SNe 
are only discovered in a fraction of the volume correspond- 
ing to the smallest distance. Furthermore, a good understand- 
ing of the detection efficiency is required to correct the small 
statistics of dim CC SNe. For example, iRichardson et al.l 
(12002) suggests that the fraction of CC SNe dimmer than 
Mb = —15 is likely larger than 20% even though their nom- 
inal luminosity function suggests a value closer to 10%, be- 
cause of the bias against dim SNe and small number statistics. 
On the other hand, we limit ourselves to a nominal distance 
of 10 Mpc and the past decade from 2000 to 2009 to study 
dim CC SNe. The majority of SN searches in this very local 
volume have limiting magnitudes of 18 mag, corresponding 
to « —12 mag. Therefore, our very local sample is only lim- 
ited by cadence and incompleteness of the galaxies observed. 
This may also explain why the dim fraction is lower at larger 
distances: by 30-40 Mpc, the limiting magnitudes of most 
searches are already more luminous than —15 mag. Finally, 
we add that even so, our results should be treated as lower 
limits, since many of the CC SNe have been discovered by 
amateurs with different systematics, galaxy samples, and ca- 
dences, i.e., not systematically collected. 

In summary, our estimate of /dim is significantly larger than 
previous studies. If possible SN impostors are not true CC 
SNe, we estimate a value for /dim that is somewhat larger 
than the SN luminosity function of LOSS. 

4.4. The Nature of Possible SN Impostors and Alternate 
Solutions to the Supernova Rate Problem? 

The dimmest objects in Table|2](SN 2002bu, SN 2008S, and 
NGC300-OT) are a special kind of explosion. Their nature — 
whether an extremely dim CC SNe or an extremely luminous 
stellar outburst — ^remains debated. Here, we discuss implica- 
tions for the supernova rate problem if they are stellar out- 
bursts. 

There is currently not enough data to conc l ude th e true na- 
ture of possible SN impostors. ISmith et akl (12009 ^ interpret 
SN 2008S as a super-Eddington wind of a ~ 2OM0 star, based 
on among o thers the lack of a ^^Co ^^Fe decay tail. On the 
other hand, iBotticella et al.l (l2009b identify in their pseudo- 
bolometric late-time light curve a decay slope consistent with 
^^Co decay, concluding an electron-capture CC SNe of a 
^ 9M0 star. Measuring the bolometric light curve is a dif- 
ficult multi-wavelength task, and further observations are re- 
quired to settle the case. The progenitor of SN 2008S has been 
identified on Spitzer pre-exp losion i mages to be a highly dust- 
enshrouded ~ IOMq star dPrieto et al...2008. .2010.) . Such 



stars are extremely rare (iKhan et al.ll2010h . suggesting it may 
be a co mmon but short-lived ph ase prior to some kind of ex- 
plosion ("Thompson et al."2009 |). The situations are similar 
for NGC30 0-OT (Bond et al..,2009t iBerger et al.l.2q09) and 
SN 2002bu (iThompson et al.ll2009HSmifli et al.ll201 ih . 

If possible SN impostors are not CC SNe, the dim frac- 
tion is only somewhat larger than the luminosity function of 
LOSS, and missing dim CC SNe provide only a partial ex- 
planation to the supernova rate problem. In this case, there 
are several interpretations of the supernova rate problem, each 
with their own tensions with other measurements. We detail 
some of them below. 

First, the supernova rate problem may be interpreted as a 
large number of undetected CC SNe. For example, the frac- 
tion ofopticallydark CC SNe could be high (see, e.g.. Figure 
3 of lLien et al.ll2010l for how the optically-dark CC SNe frac- 
tion relates to the SFR and measured SNR). Dai-k CC SNe as 
the solution to the supernova rate problem is only barely al- 
lowed by neutrino considerations dLien et alJ 2010l) and may 
be im plausible on a theoretical basis (e.g.. ^ O'Connor & Otfl 
I2OI lb . An alternate possibility is that more CC SNe than 
expected are driven below survey sensitivities by heavy host 
dust attenuation. However, this would appear difficult since 
such a change in dust properties would also increase the de- 
rived SFR and hence the predicted SNR, thus not necessar- 
ily helping solve the supe rnova rate problem. As pointed out 
by IBotticella et al.l (l2008l) . the consistency of the extinction 
model used in SFR and SNR measurements need to be fur- 
ther examined. 

A second possibility is that CC SNe in small galaxies play 
a significantly larger role than expected. As discussed in Sec- 
tion [331 the efficiency of forming CC SNe in dwarf galaxies 
may be higher than in normal galaxies. If the enhancement 
exists in a wider range of small galaxies, the effect of miss- 
ing small galaxies could be larger than we estimate. How- 
ever, there is a potential caveat with the CC SNe in small 
galaxies for helping solve the supernova rate problem. It is 
convenient to re-express the supernova rate problem in terms 
of the efficiency, or A^ccsn /M, the number of CC SNe per 
unit stellar mass formed. The measured cosmic SNR and 
the measured cosmic SFR imply O.OO41/M0, whereas the 
nominal prediction is O.OO88/M0 (see Section |2l); the fac- 
tor of ^ 2 difference is another statement of the supernova 
rate problem. If the efficiency in small galaxies is taken to 
be enhanced, the observed A^ccsn is corrected upwards, so 
that the observed iVccsN/-A^ is closer to nominal predictions. 
However, the same enhanced efficiency will also increase the 
Nqc^^/M used in predictions, unless the efficiency in nor- 
mal galaxie s is sub-nominal. This is not seen in the LOSS 
galaxies, as iMaoz etalj (1201 ih find A^ccsn/M = 0.0059- 
O.OIO/M0 for a subset of LOSS galaxies with SDSS specti'o- 
scopic SFR measurements. Correcting for systematic effects, 
their best estimate is Ncc/M = 0.010 ± 0.002/Mq, close 
to the nominal value. Therefore, enhanced CC SNe formation 
efficiencies in small galaxies may not bridge the supernova 
rate normalization discrepancy. Whether this is the case re- 
mains to be investigated with more data. For example, the 
known aperture bias from the finite SDSS fiber diameter is 
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most c ritical for low-redshif t galaxies like those in the LOSS 
survey dKewlev et alj|2005l) . This would tend to underesti- 
mate the SFR in the LOSS galaxies, potentially leading to an 
artificially larger Nqcsn / M . This may allow some parameter 
space for CC SNe in small galaxies to help bridge the cosmic 
SNR normalization mismatch. 

The final possibility is that the cosmic SFR is overesti- 
mated. In particular, the observed efficiency in the LOSS 
galaxies, Nqcsn/M = 0.010±0.002/Mo, implies that there 
is very little room for missing CC SNe in the LOSS galax- 
ies. The supernova rate problem could then be interpreted as 
an overestimated cosmic SFR. However, this may not help 
solve the supernova rate problem. Say there was some new 
effect causing the cosmic SFR to be overestimated, for exam- 
ple contamination of SFR estimators or systematic offset in 
SFR calibrations. The same effect would also cause the SFR 
in the LOSS galaxies to be overestimated, because the same 
SFR estimators are used. Thus, while reducing the cosmic 
SFR will increase the observed efficiency closer to nominal 
predictions, the efficiency in the LOSS galaxies will also in- 
crease to almost twice nominal predictions O.O2/M0). In 
other words, the supernova rate problem has only shifted. 

Additionally, the nominal uncertainties in the SFR input 
physics do not easily tolerate a factor of 2 reduction in the 
cosmic SFR (Section[TT)- The largest uncertainty in the mea- 
sured SFR is due to dust correction. However, to reduce the 
cosmic SFR sufficiently requires removing most of the dust 
correction. Apart from the need to redefine our understanding 
of the FIR universe, demanding no dust would also reduce the 
SNR measurements. The other, and perhaps more intriguing, 
possibility, is the presence of previously unknown contami- 
nations to the SFR indicators. However, several SFR indica- 
tors used must be systematically contaminated by the same 
degree, making this seem unlikely. In either case, the impli- 
cations for affecting the SFR estimators would go beyond the 
supernova rate problem to a wide range of topics. 

5. SUMMARY AND IMPLICATIONS 

We find that the SNR predicted from the cosmic SFR, us- 
ing nominal assumptions about the occurrence efficiency of 
optically luminous CC SNe, is higher than the SNR measured 
by SN surveys. Fits to the data demonstrate a normalization 
discrepancy of a factor of ^ 2 at the ^ 2a level: we term 
this the "supernova rate problem." At high redshifts {z > 0.4) 
the measurement uncertainties are large and the significance 
is lower We systematically explore resolutions. First, we in- 
vestigate the various inputs and assumptions that feed into the 
SNR predictions and show that their uncertainties do not tol- 
erate a normalization reduction by a factor of ^ 2. We then 
explore dim CC SNe (dimmer than — 15 mag) that would have 
been difficult to discover in SNR measurements. From a sam- 
ple of very local (< 10 Mpc) CC SNe, we estimate that as 
many as ^ 50% of CC SNe may be dim, which could largely 
explain the SNR discrepancy. Our analysis is enabled by the 
recent precise measurements of the cosmic SFR and the SNR. 
In particular, the extensive SN survey by LOSS greatly con- 
tributed to the discovery of local CC SNe. 

A large dim fraction has been suggested before bv lSchaefeil 



d 19961) based on historical CC SNe. Here, we analyze more 
modern data and assess the completeness of the CC SNe 
sample. We find that the dim fraction is significantly larger 
than current luminosity functions used in SFR measurements. 
However, the exact dim fraction is still uncertain for several 
reasons: (1) the nature of the dimmest CC SNe are ambigu- 
ous at present. They appear to be unlike previous classified 
explosions and, if they are not true CC SNe, the dim frac- 
tion is more modest, closer to ~ 30%; (2) dim CC SNe have 
not been systematically searched for, so many may have been 
missed. Until a systematic survey addresses this, our esti- 
mates for dim CC SNe should be treated as lower limits; (3) 
based on the current statistics, the dim fraction follows a fairly 
broad probability distribution; (4) finally, CC SNe could be 
missed for other reasons. For example, a fraction of CC SNe 
associated with the prompt formation of black holes could be 
even dimmer and not observable by current optical, IR, or ra- 
dio surveys. 

The supernova rate problem raises some interesting impli- 
cations. We are left with three main possible outcomes. 

• First, the fraction of dim CC SNe is high (/>-i5 
30%-50%). This is likely if the possible SN impostors 
are true CC SNe. The supernova rate problem is solved, 
with the implication that half of stars with masses 8- 
40 Mq are producing dim CC SNe, either due to dust 
obscuration or being intrinsically weak. 

• Second, if possible SN impostors are not CC SNe, the 
dim fraction is />-i5 ^ 20-40%, i.e., only slightly 
higher than the most re cent SN luminosity function of 
LOSS jLi et al.ll20l"Tbl) . The majority of dim CC SNe 
are due to heavily obscured CC SNe. The supernova 
rate problem could be explained by a high fraction of 
optically dark CC SNe, although the required fraction 
is somewhat higher th an theoretical expectations (e.g., 
IO'Connor&Ottll20TTl) . 

• Third, the dim fraction is />_i5 20%-40%, and 
the supernova rate problem is explained by systematic 
changes in our understanding of star formation or CC 
SN formation. For example, the efficiency of forming 
CC SNe in small galaxies may be higher than in normal 
galaxies. However, this may not help solve the super- 
nova rate problem because the predicted SNR is also 
increased. Another example is that SFR estimators may 
be systematically overestimating the true cosmic SFR, 
although any explanation must go beyond current un- 
derstanding of SFR estimation. If this turns out to be 
the cause, the implications would go beyond the super- 
nova rate problem to a wide range of topics. 

To resolve these outcomes, detailed studies of the very local 
(< 10 Mpc) SFR using various indicators are needed to set the 
expected normalization of all CC SNe including optically lu- 
minous, dim, and dark CC SNe. The increased coverage of the 
Southern hemisphere by e.g., the CHilean Automatic Super- 
nova sEarch (CHASE; 'Pi gnata et al.l2009l ) will aid by provid- 
ing a more complete discovery list of SNe. To verify the dim 
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fraction, we need to capture all collapsing massive stars: op- 
tically luminous, dim, and dark. The All-Sky Aut omated Sur- 
vey for the B rightest Supernovae (ASAS-SN; e.g., 'Khan et alj 
bOllHStoUet al. 201 1) will find CC SNe in a volume-limited 
sample of nearby galaxies (< 30 Mpc) that will contribute to 
a better census of relatively low-luminosity events. In fact, 
one previously missed CC SN has already been identified in 
ASAS and CRTS images jPrieto et al.ll201lh . The PTF will 
collect a large number of early-discovered CC SNe and allow 
a hig h-statistics study of CC SNe varieties and host prope r- 



ties dLaw et al.l2009tlLien & Fieldsl2009l:lArcavi et al.l 2010). 
Already, dim CC SNe (potentially dimmer than SN 1999br 
and PTFlOvdl/S N 2010dl) have been discovered by the PTF 
-Yaml20TTI) . Another interesting possibility is to observe 
the disappearance of massive stars as opposed to the appear- 
ance of_timrsxj^\o^ons^hich will reveal the dark collapse 
rate ("Koch anek et al.ll2008f ). Progenitor studies will be impor- 
tant in connecting future SFR and CC SNe studies. 

It will ultimately become possible to probe all kinds of 
collapsing stars with the deployment of next-generation non- 
electro magnetic probes su ch as neutrinos and gravitational 
waves ( iArnaud et al.ll2004 . All core collapses — whether op- 
tically luminous or dim or dark — make neutrinos in compara- 
ble numbers. Proposed megaton neutrino detectors would be 
sensitive to core collapses in the several Mp c volume, reveal- 
ing the true collap se rate of massive stars ( Ando et al1l2005[ 
iKistler et alj l2008l) . The nature of possible SN impostors 
could also be tested with the observation or non-observation 
of a neutrino signal. 

The diffuse SN neutrino back ground (DSNB; see, e.g., 
iBeacomI I2OI0I: iLunardinil l2010"l similarly probes all core 
collapses and will test the true cosmic collapse rate 
of mass ive stars. The S uper-Kamiokande limit on the 
DSNB jMalek et alJ l2003h is very near theoretical pre- 
dictions using a co smic SFR similar t o the one as- 
sumed in this paper (iHoriuchi et al.l l2009h . and prospects 



for detecting the DSNB will be significantly increased 
by the near-future Gadolinium-enhanced Super-Kamiokande 
dBeacom & VaginsI l2004t IWatanabj|200 8). We emphasize 
that while the luminous SNR in cosmic surveys is lower than 
expected, there is no evidence for a reduction in the pre- 
dicted DSNB flux, which depends on all core collapses. If the 
neutrino emission from a collapse to a black hole is signifi- 
cantly higher in energy compared t o the canonical colla pse to 
a neutron star (e.g., as discussed in lNakazato et al.|[2008i) . the 
DSNB will be sensitive to the fraction of black hole form ing 
collapses (e.g.. lLien et alJl2O10l : lKeehn & Lunardinill2010h . 

Using the full array of SN surveys, progenitor studies, sur- 
vey of disappearing stars, DSNB searches, and upcoming 
gravitational wave detectors, a comprehensive understanding 
of the deaths of massive stars and their associated observ- 
able transient phenomena will come within reach in the com- 
ing years. This will be an important complement to stud- 
ies of the progenitors of SN la and their explosion mecha- 
nisms that will be opened by similar optical SN studies to- 
gether with next-generation MeV gamma-ray detectors (e.g., 
Horiuchi & Beacom 2010). 

We thank Berto Monard for sharing photometry of 
SN 2008eh. We thank Matt Bothwell, Tim Eifler, Alex Fil- 
ippenko, Avishay Gal- Yam, Robert Kennicutt, Matt Kistler, 
Weidong Li, Amy Lien, Dan Maoz, Michael Mortenson, Jim 
Rich, and Stephen Smartt for discussions. This research made 
use of the lAU Central Bureau for Astronomical Telegrams 
and the Sternberg Astronomical Institute supernova catalogs 
and the NASA/IPAC Extragalactic Database (NED), which is 
operated by JPL/Caltech, under contract with NASA. SH and 
JFB were supported by NSF CAREER Grant PHY-0547102 
(to JFB); CSK, KZS, and TAT by NSF Grant AST-0908816; 
and JLP by NASA through Hubble Fellowship Grant HF- 
5 126 1.01 -A awarded by STScl, which is operated by AURA, 
Inc. for NASA, under contract NAS 5-2655. 



REFERENCES 



Abdo, A. A., et al. 2010, ApJ, 723, 1082 
Aharonian, E, et al. 2006, Nature, 440, 1018 
Albert, J., et al. 2008, Science, 320, 5884 
Anderson, J. R, & James, R A. 2008, MNRAS, 390, 1527 
Anderson, J. R, & James, R A. 2009, MNRAS, 399, 559 
Ando, S., Beacom, J. E, Yiiksel, H. 2005, Physical Review Letters, 95, 
171101 

Arbour, R., & Boles, T. 2003, lAU Circ, 8205, 1 

Arcavi, I., et al. 2010, ApJ, 721, 777 

Amaud, N., et al. 2004, Astroparticle Physics, 21, 201 

Amett, W. D. 1987, ApJ, 319, 136 

Amett, W. D. 1980, ApJ, 237, 541 

Amouts, S., et al. 2005, ApJ, 619, L43 

Avni, Y. 1976, ApJ, 210, 642 

Baldry, I. K., & Glazebrook, K. 2003, ApJ, 593, 258 
Baldry, I. K., et al. 2005, MNRAS, 358, 441 

Bartunov, O. S., Tsvetkov, D. ¥., & Pavlyuk, N. N. 2007, HighUghts of 

Astronomy, 14, 316 
Bazin, G., et al. 2009, A&A, 499, 653 

Beacom, J. E, & Vagins, M. R. 2004, Physical Review Letters, 93, 171101 
Beacom, J. E. 2010, Annual Review of Nuclear and Particle Science, 60, 439 
Bell, E. E 2003, ApJ, 586, 794 

Bembrick, C, Pearce, A., & Evans, R. 2002, lAU Circ, 7804, 2 
Berger, E., et al. 2009, ApJ, 699, 1850 

Biggs, J., Monard, L. A. G., & Africa, S. 2003, lAU Circ, 8236, 3 



Bond, H. E., Bedin, L. R., Bonanos, A. Z., Humphreys, R. M., Monard, 

L. A. G. B., Prieto, J. L., & Walter, E M. 2009, ApJ, 695, L154 
Bothwell, M. S., et al. 2011, arXiv: 1104.0929 
Botticella, M. T., et al. 2008, A&A, 479, 49 
Botticella, M. T., et al. 2009, MNRAS, 398, 1041 

Brinchmann, J., Chariot, S., White, S. D. M., Tremonti, C, Kauffmann, G., 

Heckman, T., & Brinkmann, J. 2004, MNRAS, 351, 1151 
Branthaler, A., et al. 2010, A&A, 516, A27 
Cappellaro, E., Evans, R., & Turatto, M. 1999, A&A, 351, 459 
Cappellaro, E., et al. 2005, A&A, 430, 83 

Cardelh, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Condon, J. J., Cotton, W. D., & Broderick, J. J. 2002, AJ, 124, 675 
Dahlen, T., et al. 2004, ApJ, 613, 189 

Dahlen, T., Strolger, L., & Riess, A. G. 2010, Bulletin of the American 

Astronomical Society, 42, 360 
Delia Valle, M., et al. 2006, Nature, 444, 1050 
Elias-Rosa, N., et al. 2010, ApJ, 714, L254 
Essey, W., & Kusenko, A. 2010, Astroparticle Physics, 33, 81 
Essey, W., Kalashev, O., Kusenko, A., & Beacom, J. E 2011, ApJ, 731, 51 
Eardal, M. A., Katz, N., Weinberg, D. H., & Dave, R. 2007, MNRAS, 379, 

985 

Eerrara, A., & Ricotti, M. 2006, MNRAS, 373, 571 
Eioc, M., & Rocca-Volmerange, B. 1997, A&A, 326, 950 
Eraser, M., et al. 2010a, ApJ, 714, L280 
Eraser, M., et al. 2010b, arXiv: 101 1.6558 
Eruchter, A. S., et al. 2006, Nature, 441, 463 



The supernova rate does not match the star formation rate 



17 



Fiyer, C. L. 1999, ApJ, 522, 413 

Fynbo, J. P. U., et al. 2006, Nature, 444, 1047 

Gal- Yam, A., Ofek, E. O., & Shemmer, O. 2002, MNRAS, 332, L73 

Gal- Yam, A., et al. 2007, ApJ, 656, 372 

Gal- Yam, A., & Leonai'd, D. C. 2009, Nature, 458, 865 

Gal- Yam, A., et al. 201 1, ApJ, submitted 

Gal- Yam, A. 20 11, private communication 

Gallego, J., Zamorano, J., Aragon-Salamanca, A., & Rego, M. 1995, ApJ, 
455, LI 

Giavalisco, M., et al. 2004, ApJ, 600, L103 

Graur, O., et al. 2011, arXiv:l 102.0005 

Hamuy, M. 2003, ApJ, 582, 905 

Hanish, D. J., et al. 2006, ApJ, 649, 150 

Hatano, K., Branch, D., & Deaton, J. 1998, ApJ, 502, 177 

Hauser, M. G., & Dwek, E. 2001, ARA&A, 39, 249 

Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Hartmann, D. H. 

2003, ApJ, 591,288 
Hopkins, A. M. 2004, ApJ, 615, 209 
Hopkins, A. M., & Beacom, J. E 2006, ApJ, 651, 142 
Hopkins, P. E, Hemquist, L., Cox, T. J., Di Matteo, T., Robertson, B., & 

Spiingel, V. 2006, ApJS, 163, 1 
Horiuchi, S., Beacom, J. F, & Dwek, E. 2009, Phys. Rev. D, 79, 083013 
Horiuchi, S., & Beacom, J. F 2010, ApJ, 723, 329 
Hornoch, K. 2002, lAU Circ, 7923, 5 
Iglesias-Paramo, J., et al. 2006, ApJS, 164, 38 
Jacques, C, & Pimentel, E. 2005, lAU Circ, 8482, 1 
James, R A., Knapen, J. H., Shane, N. S., Baldry, I. K., & de Jong, R. S. 

2008, A&A, 482, 507 
Janka, H.-T., Langanke, K., Marek, A., Marti'nez-Pinedo, G., Muumlller, B. 

2007, Phys. Rep., 442, 38 
Jarosik, N., et al. 2011, ApJS, 192, 14 

Kalirai, J. S., Hansen, B. M. S., Kelson, D. D., Reitzel, D. B., Rich, R. M., & 

Richer, H. B. 2008, ApJ, 676, 594 
Kankare, E., et al. 2008, ApJ, 689, L97 

Karachentsev, I. D., Karachentseva, V. E., Huchtmeier, W. K., & Makarov, 

D. 1.2004, AJ, 127, 2031 

Keehn, J. G., & Lunardini, C. 2010, arXiv: 1 012.12741 
Kennicutt, R. C, Jr. 1998, ARA&A, 36, 189 

Kennicutt, R. C, Jr., Lee, J. C, Funes, S. J., Jose G., Sakai, S., & Akiyama, 

S. 2008, ApJS, 178, 247 
Kewley, L. J., Geller, M. J., Jansen, R. A., & Dopita, M. A. 2002, AJ, 124, 

3135 

Kewley, L. J., Jansen, R. A., & Geller, M. J. 2005, PASP, 117, 227 

Khan, R., Stanek, K. Z., Prieto, J. L., Kochanek, C. S., Thompson, T. A., & 

Beacom, J. F 2010, ApJ, 715, 1094 
Khan, R., et al. 201 1, ApJ, 726, 106 

Kisder, M. D., Yuksel, H., Ando, S., Beacom, J. F, & Suzuki, Y. 2008, 
arXiv:0810.1959 

Kochanek, C. S., Beacom, J. F, Kisder, M. D., Prieto, J. L., Stanek, K. Z., 

Thompson, T. A., Yuksel, H. 2008, ApJ, 684, 1336 
Kochanek, C. S., Szczygiel, D. M., & Stanek, K. Z. 2010, arXiv: 1010.3704 
Kotake, K., Sato, K., & Takahashi, K. 2006, Reports on Progress in Physics, 

69, 971 

Kulkarni, S. R., et al. 2007, Nature, 447, 458 
Lagache, G., Puget, J.-L., & Dole, H. 2005, ARA&A, 43, 727 
Lattimer, J. M., & Prakash, M. 2007, Phys. Rep., 442, 109 
Law, N. M., et al. 2009, PASP, 121, 1395 

Le Delhou, M., Lacey, C. G., Baugh, C. M., & Morris, S. L. 2006, MNRAS, 
365, 712 

Le Floc'h, E., et al. 2003, A&A, 400, 499 
Le Floc'h, E., et al. 2005, ApJ, 632, 169 

Leaman, J., Li, W., Chomock, R., & Filippenko, A. V. 2011, MNRAS, 412, 
1419 

Li, W., et al. 2011, MNRAS, 412, 1441 

Li, W., Chomock, R., Leaman, J., Filippenko, A. V., Poznanski, D., Wang, 
X., Ganeshalingam, M., & Mannucci, F. 2011, MNRAS, 412, 1473 

Ly, C, Lee, J. C, Dale, D. A., Momcheva, I., Salim, S., Staudaher, S., 
Moore, C. A., & Finn, R. 2011, ApJ, 726, 109 

Lien, A., & Fields, B. D. 2009, JCAR 1, 47 

Lien, A., Fields, B. D., & Beacom, J. F. 2010, Phys. Rev. D, 81, 083001 
Lonsdale Persson, C. J., & Helou, G. 1987, ApJ, 314, 513 
Lunardini, C. 2010, arXiv: 1007.3252 

Madau, P, Pozzetd, L., & Dickinson, M. 1998, ApJ, 498, 106 
Malek, M., et al. 2003, Physical Review Letters, 90, 061 101 
Mannucci, E, et al. 2003, A&A, 401, 519 

Mannucci, E, Delia Valle, M., & Panagia, N. 2007, MNRAS, 377, 1229 



Maoz, D., Mannucci, F, Li, W., Filippenko, A. V., Delia Valle, M., & 

Panagia, N. 2011, MNRAS, 412, 1508 
Martin, D. C, et al. 2005, ApJ, 619, L59 

Martin, R., Yamaoka, H., Monard, L. A. G., & Africa, S. 2005, lAU Circ, 
8496, 1 

Mattila, S., Meikle, W. P. S., Groeningsson, P., Greimel, R., Schinner, M., 

Acosta-Pulido, J. A., & Li, W. 2004, lAU Circ, 8299, 2 
Massey, P, Waterhouse, E., & DeGioia-Eastwood, K. 2000, AJ, 1 19, 2214 
Massey, P, DeGioia-Eastwood, K., & Waterhouse, E. 2001, AJ, 121, 1050 
Matteucci, F, & Greggio, L. 1986, A&A, 154, 279 
Maund, J. R., et al. 2006, MNRAS, 369, 390 
McNaught, R. H. 2003, lAU Circ, 8152, 3 

Miyaji, S., Nomoto, K., Yokoi, K., & Sugimoto, D. 1980, PASJ, 32, 303 
Monard, L. A. G. 2008a, Central Bureau Electronic Telegrams, 1315, 1 
Monard, L. A. G. 2008b, Central Bureau Electronic Telegrams, 1445, 1 
Monard, L. A. G. 2009, Central Bureau Electronic Telegrams, 1867, 1 
Morrell, N., & Stritzinger, M. 2008, Central Bureau Electronic Telegrams, 
1335, 1 

Nagamine, K., Cen, R., Hernquist, L., Ostiiker, J. P., & Springel, V. 2004, 
ApJ, 610, 45 

Nagashima, M., Lacey, C. G., Okamoto, T, Baugh, C. M., Erenk, C. S., & 

Cole, S. 2005, MNRAS, 363, L31 
Nakazato, K., Sumiyoshi, K., Suzuki, H., & Yamada, S. 2008, Phys. Rev. D, 

78, 083014 
Nomoto, K. 1984, ApJ, 277, 791 
O'Connor, E., & Ott, C. D. 201 1, ApJ, 730, 70 
OiT, M. R., Krennrich, F, & Dwek, E. 201 1, ApJ, 733, 77 
Pascale, E., et al. 2009, ApJ, 707, 1740 
Pastorello, A., et al. 2004, MNRAS, 347, 74 
Pastorello, A., et al. 2006, MNRAS, 370, 1752 
Pastorello, A., et al. 2007, Nature, 449, 1 
Pastorello, A., et al. 2008, MNRAS, 389, 955 

Patat, E, Barbon, R., Cappellaro, E., & Turatto, M. 1994, A&A, 282, 731 
Peacock, J. A., & Dodds, S. J. 1994, MNRAS, 267, 1020 
Perez-Gonzalez, P. G., Zamorano, J., Gallego, J., Aragon-Salamanca, A., & 

Gil de Paz, A. 2003, ApJ, 591, 827 
Pignata, G. 2009, Central Bureau Electronic Telegrams, 1902, 1 
Pignata, G., et al. 2009, American Institute of Physics Conference Series, 

1111, 551 
Prieto, J. L., et al. 2008, ApJ, 681, L9 

Prieto, J. L., Szczygiel, D. M., Kochanek, C. S., Stanek, K. Z., Thompson, 
T. A., Beacom, J. E, Garnavich, P M., & Woodward, C. E. 2010, 
arXiv: 1007.0011 

Prieto, J. L., et al. 201 1, in prep. 

Poelarends, A. J. T, Herwig, E, Langer, N., & Heger, A. 2008, ApJ, 675, 
614 

Pojmanski, G. 2007, lAU Circ, 8875, 1 
Pozzo, M., et al. 2006, MNRAS, 368, 1169 
Raffelt, G. G. 1990, Phys. Rep., 198, 1 
Raffelt, G. G. 2000, Phys. Rep., 333, 593 
Rich, J. 2011, private communication 

Richardson, D., Branch, D., Casebeer, D., Millard, J., Thomas, R. C, & 

Baron, E. 2002, AJ, 123, 745 
Rujopakam, W., et al. 2010, ApJ, 718, 1171 

Sahu, D. K., Anupama, G. C, Srividya, S., & Muneer, S. 2006, MNRAS, 

372, 1315 
Sahm, S., et al. 2007, ApJS, 173, 267 
Schaefer, B. E. 1996, ApJ, 464, 404 
Schiminovich, D., et al. 2005, ApJ, 619, L47 
Schiminovich, D., et al. 2007, ApJS, 173, 315 
Schlegel, D. J., Einkbeiner, D. R, & Davis, M. 1998, ApJ, 500, 525 
Serjeant, S., Grappioni, C, & Oliver, S. 2002, MNRAS, 330, 621 
Silva, L., Maiolino, R., & Granato, G. L. 2004, MNRAS, 355, 973 
Singer, D., Pugh, H., & Li, W. 2004, lAU Circ, 8297, 2 
Smartt, S. J. 2009, ARA&A, 47, 63 

Smartt, S. J., Eldridge, J. J., Crockett, R. M., & Maund, J. R. 2009, MNRAS, 

395, 1409 
Smith, N., et al. 2009, ApJ, 697, L49 

Smith, N., Li, W., Silvemian, J. M., Ganeshalingam, M., & Filippenko, A. V. 

2010, arXiv:1010.3718 
Smith, N., et al. 201 1, ApJ, 732, 63 
Soifer, B. T, & Neugebauer, G. 1991, AJ, 101, 354 

Somerville, R. S., Primack, J. R., & Faber, S. M. 2001, MNRAS, 320, 504 
StoU, R., Prieto, J. L., Stanek, K. Z., Pogge, R. W., Szczygiel, D. M., 

Pojmanski, G., Antognini, J., & Yan, H. 201 1, ApJ, 730, 34 
Taylor, M., Cinabro, D., & SDSS-Il Supernova Survey Team 2009, Bulletin 

of the American Astronomical Society, 41, 252 



18 



Horiuchi et al. 



Thompson, T. A., Prieto, J. L., Stanek, K. Z., Kistler, M. D., Beacom, J. F., 

& Kochanek, C. S. 2009, ApJ, 705, 1364 
Todini, P., & Ferrara, A. 2001, MNRAS, 325, 726 
Tully, R. B., Rizzi, L., Shaya, E. J., Courtois, H. M., Makarov, D. I., & 

Jacobs, B. A. 2009, AJ, 138, 323 
Turatto, M., et al. 1998, ApJ, 498, L129 
Valenti, S., et al. 2008, ApJ, 673, L155 

Van Dyk, S. D., Filippenko, A. V., Chomock, R., Li, W., & Challis, R M. 

2005, RASP, 117, 553 
Vink6, J., et al. 2006, MNRAS, 369, 1780 
Watanabe, H. 2008, International Cosmic Ray Conference, 5, 1421 
Watson, C. R., et al. 2009, ApJ, 696, 2206 

Wilkins, S. M., Hopkins, A. M., Trentham, N., & Tojeiro, R. 2008, MNRAS, 
391, 363 

Wilkins, S. M., Trentham, N., & Hopkins, A. M. 2008, MNRAS, 385, 687 



Williams, K. A., Bolte, M., & Koester, D. 2009, ApJ, 693, 355 
Wolf, C, Meisenheimer, K., Rix, H.-W, Borch, A., Dye, S., & 

Kleinheinrich, M. 2003, A&A, 401, 73 
Woosley, S. E. 1993, ApJ, 405, 273 

Young, D. R., Smartt, S. J., Manila, S., Tanvir, N. R., Bersier, D., Chambers, 

K. C, Kaiser, N., & Tonry, J. L. 2008, A&A, 489, 359 
Yiiksel, H., Kistler, M. D., Beacom, J. F., & Hopkins, A. M. 2008, ApJ, 683, 

L5 

Yun, M. S., Reddy, N. A., & Condon, J. J. 2001, ApJ, 554, 803 
Zampieri, L., PastoreUo, A., Turatto, M., CappeUaro, E., Benetti, S., 

Altavffla, G., MazzaH, R, & Hamuy, M. 2003, MNRAS, 338, 711 
Zhang, T., Wang, X., Li, W, Zhou, X., Ma, J., Jiang, Z., & Chen, J. 2006, 

AJ, 131, 2245 



